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I.  INTRODUCTION 

This  is  the  annual  report  and  summary  of  accomplishments  for  the  second  year  of 
a  three-yeax  program  of  research  on  laser  physics  and  laser  techniques  carried  out  in 
the  research  group  of  Professor  A.  E.  Siegman  at  Stanford  University  and  supported 
by  AFOSR  Contract  No.  F49620-S9-K-0004,  covering  the  period  from  15  March  1990 
to  15  March  1991. 

This  report  briefly  summarizes  results  on  laser  resonators  and  laser  resonator 
and  optical  beam  software;  new  techniques  for  laser  beam  characterization  and  laser 
beam  quality  measurement;  the  generation  of  ultrashort  optical  pulses  in  the  infrared; 
and  the  meaisurement  of  a  fundamental  excess  quantum  noise  effect  in  ultrastable 
laser  oscillators.  Accomplishments  in  each  of  the  areas  of  research  covered  by  this 
program  are  described  in  more  detail  in  the  following  sections,  and  in  various  reprints 
or  preprints  attached  as  Appendices 


II.  RAY-PULSE  MATRIX  ANALYSIS 

A  very  useful  ray-pulse  matrix  approach  to  optical  pulse  propagation  was  devel¬ 
oped  earlier  under  this  program  by  Dr.  A.  Kostenbauder,  as  described  in  last  year’s 
annual  report.  This  analysis  appeared  in  print  during  the  current  period  in  the  journal 
article; 

A.  G.  Kostenbauder,  “Ray -pulse  matrices;  A  rational  treatment  for 
dispersive  optical  systems,”  IEEE  J.  Quantum  Electron.  QE— 26. 
1148-1157  (June  1990). 

This  particular  project  is  now  completed;  and  a  reprint  of  this  publication  is  attached 
as  Appendix  A. 
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III.  OPTICAL  BEAM  AND  RESONATOR  CALCULATIONS 

Our  group  has  made  many  contributions  in  previous  years  to  optical  resonator 
concepts,  and  has  also  done  much  work  on  fundamental  methods  for  the  analysis  of 
optical  beam  propagation  and  the  calculation  of  laser  resonator  modes;  and  these 
topics  continue  to  be  a  productive  area  of  work  in  the  group. 

In  particular  during  the  current  reporting  period  we  have  completed  the  devel¬ 
opment  of  an  integrated  package  of  software  programs  with  which  optics  engineers 
and  laser  researchers  can  carry  out  optical  beam  propagation  and  resonator  mode  cal¬ 
culations  in  real  time  on  personal  desktop  computers.  This  package,  which  we  have 
named  PARAXIA  because  of  its  basic  dependence  on  the  paraxial  wave  approximation, 
includes: 

o  An  ABCD  program  for  carrying  out  complex-valued  ABCD  matrix  calculations, 
including  two-transverse-dimension  or  astigmatic  systems.  This  program  is  par¬ 
ticularly  useful  in  both  stable  and  unstable  optical  resonator  evaluation  and  in 
optical  system  design. 

o  A  program  FRESNEL  which  carries  out  accurate  optical  beam  propagation  and 
diffraction  calculations  in  rectangular  or  radial  coordinates  using  fast  transform 
algorithms.  This  program  also  includes  scripting,  iteration,  and  external  function 
capabilities  for  doing  iterative  Fox-and-Li  resonator  mode  calculations. 

o  A  program  VSOURCE  which  implements  the  fast  and  efficient  virtual  source 
algorithm  for  calculating  the  modes  of  rectangular  and  circular  unstable  laser 
resonators,  including  both  lowest  and  higher-order  modes. 

o  A  VRM  program  for  designing  and  evaluating  the  mode  properties  of  gaussian 
variable- reflectivity-mirror  resonators. 

This  software  package  was  initially  made  available  to  other  university,  industrial  and 
government  laser  laboratories  through  the  Stanford  Software  Distribution  Center,  and 
at  least  several  dozen  copies  were  distributed  to  other  laser  researchers  in  this  fashion, 
including  site  licenses  to  several  large  laser  Arms.  A  brochure  describing  this  pack¬ 
age  and  its  capabilities  is  attached  as  Appendix  B.  We  have  received  very  favorable 
comments  concerning  the  usefulness  of  this  software  from  a  number  of  university  and 
industrial  users.  The  package  has  now  been  licensed  to  Genesee  Optics  Software.  Inc., 
of  Rochester.  NY,  a  widely  known  vendor  of  lens  design  and  optics  software  packages 
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for  the  Macintosh,  with  the  expectation  that  it  will  receive  still  wider  distribution  and 
use  in  the  future  through  this  channel. 

These  programs  were  initially  written  primarily  for  our  internal  use  on  this  project, 
and  we  have  made  extensive  use  of  this  software  ourselves  in  our  own  resonator  design 
and  laser  beam  quality  research,  as  described  in  later  sections  of  this  report.  Two 
such  projects  which  were  carried  out  using  this  software  and  which  are  now  essentially 
complete,  although  the  results  are  still  being  prepared  for  submission,  are; 

K.  Yasui,  P.  Mussche,  J-L  Doumont,  and  A.  E.  Siegman,  ‘‘Off-axis 
one-sided  positive  and  negative  branch  unstable  resonators,”  Appl. 

Optics,  in  preparation  (1991). 

K.  Yasui,  P.  Mussche,  J-L  Doumont,  and  A.  E.  Siegman,  “Compar¬ 
ison  of  misalignment  sensitivities  for  positive  and  negative  branch 
strip  confocal  unstable  resonators,”  AppJ.  Optics,  in  preparation 
(1991). 

A  PhD  dissertation  reporting  on  the  development  of  these  programs  and  some  of 
their  applications  was  also  essentially  completed  during  the  current  period  and  will  be 
submitted  early  in  the  coming  period: 

Jean-luc  Doumont.  Laser  Beam  and  Resonator  Calculations  on 
Desktop  Computers,  Ph.D.  Dissertation,  Department  of  Applied 
Physics,  Stanford  University,  to  be  submitted  (June  1991). 

The  abstract  of  this  forthcoming  dissertation  is  attached  as  Appendix  C. 

As  a  secondary  effort  in  resonator  theory,  we  also  carried  out  during  this  reporting 
period  a  slightly  extended  analysis  and  evaluation  of  the  so-called  self-filtering  unstable 
resonator  (SFUR)  design  concept,  leading  to  a  pending  journal  article: 

A.  E.  Siegman,  “Performance  limitations  of  the  self-filtering  unsta¬ 
ble  resonator,”  Optics  Commun.,  accepted  for  publication  (1991). 

A  preprint  of  this  communication  is  attached  as  Appendix  D.  The  primary  point  to 
this  brief  communication  is  that  the  SFUR  concept,  which  has  been  widely  discussed 
by  others,  does  offer  very  good  mode  properties  for  laser  systems  with  low  Fresnel 
numbers;  but  it  is  not  a  useful  solution  for  the  more  difficult  (and  generally  more 
realistic)  czise  of  large- Fresnel- number  laser  resonators.  The  reason  is  that  as  the 
Fresnel  number  increases  much  above  a  few  times  unity,  the  output  coupling  of  the 
SFUR  design  inherently  must  become  very  large,  or  the  effective  reflectivity  of  the 
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output  mirror  becomes  very  small,  and  so  this  type  of  resonator  cannot  be  used  with 
any  normal  gain  medium. 

An  overall  survey  of  advances  in  laser  resonators  and  resonator  theory  over  the 
past  two  decades,  including  a  review  of  current  and  expected  future  developments 
in  this  field,  was  also  prepared  and  presented  as  an  invited  lecture  at  the  1990 
CLEO/QELS  meeting: 

A.  E.  Siegman,  ‘‘Laser  Resonators;  Retrospective  and  Prospective’’, 

R.V.  Pole  Memorial  Lecture.  CLEO/QELS  ’90,  Anaheim,  CA  (21 
May  1990). 

An  outline  of  this  lecture  and  associated  transparencies  is  attached  as  Appendix  E  to 
this  report. 

IV.  LASER  BEAM  QUALITY  DEFINITION  AND  MEASUREMENT 

A  major  new  area  of  research  in  our  group  during  the  present  and  the  immediately 
preceding  reporting  periods  has  been  the  general  topic  of  laser  beam  quality  definition 
and  measurement.  This  topic  is  a  major  new  effort  in  our  group,  although  it  grows 
directly  ont  of  our  etirlier  research  interests  in  optical  beams  and  resonators;  and  we 
believe  it  may  have  a  substantial  impact  on  the  development  and  improvement  of  many 
types  of  lasers. 

By  way  of  background,  since  the  earliest  days  of  the  laser  there  has  been  great 
interest  in  the  transverse  mode  properties  of  lasers,  and  much  effort  has  been  devoted 
to  obtaining  laser  output  beams  that  were  ‘‘diffraction  limited”  in  character,  since  this 
characteristic  of  a  laser  beam  has  a  critical  impact  on  the  usefulness  of  a  laser  beam, 
e.g.,  on  its  propagation  and  beam-spread  properties  for  laser  radars,  or  on  the  focusing 
properties  of  a  laser  beam  for  laser  materials  processing.  Despite  the  importance  of 
laser  beam  quality,  however,  there  has  been  very  little  work  in  the  past  either  on 
rigorous  and  meaningful  definitions  of  laser  beam  quality,  or  on  practical  tools  for 
accurate  measurement  of  these  laser  beam  characteristics. 

In  work  under  this  program  during  the  preceding  reporting  period  our  group  first 
developed  a  very  meaningful  and  widely  useful,  but  also  precise  and  universal  definition 
for  laser  beam  quality.  This  work  produced  a  rigorous  definition  for  the  so-called 
factor”  of  a  laser  beam,  based  on  the  spatial  and  angular  second  moments  or  standard 
deviations  of  a  laser  beam  and  their  space- beamwidth  product.  Our  analysis,  which 
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was  based  on  the  so-called  "moments  method"  for  laser  beam  propagation,  in  fact  gave 
a  complete  description  of  the  propagation  properties  of  any  cirbitrary  real,  distorted, 
multimode,  or  astigmatic  laser  beam  in  terms  of  only  six  easily  measurable  parameters: 
the  rezd-beam  waist  size,  the  real-beam  waist  location,  and  the  real-beam  quality  factor 
in  each  transverse  coordinate.  The  motivation  behind  this  work  was  to  provide 
quantitative,  measureable,  and  widely  useful  measures  both  of  “beam  quality"  and  of 
beam  propagation  parameters  in  the  form  of  a  few  numbers  which  can  be  attached  to 
a  given  optical  beam  and  be  valid  at  any  point  along  the  beam,  rather  than  merely 
stating  some  quantitative  measure  of  performance  such  as  far-field  brightness  in  watts 
per  steradian  or  focused  spot  size  in  microns  for  a  specific  laser  device  under  specific 
operating  conditions. 

This  earlier  analysis  focused  on  the  axial  variation  of  the  intensity  profile  of  an 
arbitrary  real  beam.  As  one  extension  of  this  analysis,  during  the  current  period  we 
were  able  to  extend  this  analysis  to  define  an  effective  spherical  radius  of  curvature 
for  an  arbitrary  laser  beam,  even  one  with  an  arbitrarily  distorted  or  wrinkled  phase 
front.  This  result  hcis  now  been  published  as; 

A.  E.  Siegman.  “Defining  the  effective  radius  of  curvature  for  a 
nonideal  optical  beam.  ’  IEEE  J.  Quantum  Electron.  QE-27,  1146- 
1148  (May  1991). 

This  new  result  has  proven  to  be  extremely  useful  in  evaluating  the  effects  on  beam 
quality  of,  for  example,  phasefront  distortion  by  spherical  aberration  or  other  phase 
aberration  effects.  A  reprint  of  this  publication  is  attached  as  Appendix  F. 

In  parallel  with  and  partly  based  on  these  theoretical  developments,  at  least  one 
commercial  firm  (Coherent,  Inc.)  has  developed  over  the  past  two  years  a  compar¬ 
atively  simple  and  practical  instrument  for  measuring  laser  beam  quality  and  beam 
propagation  parameters  in  real  time.  For  simplicity  this  particular  instrument  mea- 
smes  beam  diameter  in  terms  of  a  knife-edge  clip  width,  rather  than  a  direct  evaluation 
of  the  second  moment.  The  approximations  involved  in  this  design  compromise  are 
evaluated  in  a  joint  publication  from  our  two  laboratories: 

A.  E.  Siegman.  M.  \V.  Sasnett.  and  T.  F.  Johnston.  Jr..  “Choice 
of  clip  levels  for  beam  width  measurements  using  knife-edge  tech¬ 
niques".  IEEE  J.  Quantum  Electron.  QE-27.  1098-1104  (April 
1991). 
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a  reprint  of  which  is  attached  as  Appendix  G.  Our  group  has  also  been  using  one  of 
the  Coherent  instruments  for  various  laser  beam  evaluations  as  described  below. 

In  addition,  since  the  Coherent  instrument  operates  only  with  cw  beams,  while 
the  beam  parameters  of  pulsed  lasers  are  also  of  substantial  interest,  during  the  past 
period  we  tissembled  and  tested  in  our  laboratory  a  single-shot  pulsed  beam  quality 
and  beam  propagation  measurement  system  using  CCD-camera  arrays  and  computer 
frame- grabber  boards.  This  approach  permits  detailed  evaluation  of  the  measured 
beam  profiles  and  evaluation  of  the  beam  quality  using  various  algorithms  in  software. 
This  work  is  described  in: 

J.  A.  Ruff  and  A.  E.  Siegman,  "Single  pulse  laser  beam  quality 
measurements  using  a  CCD  camera  system,”  Appl.  Optics,  accepted 
for  publication  (1991). 

preprint  of  this  publication  is  attached  as  Appendix  H.  An  improved  version  of 
this  approach  is  expected  to  be  \ery  useful  for  several  types  of  pulsed  lasers,  such  as 
Q-switched  solid-state  lasers  and  pulsed  excimer  lasers. 

As  one  particularly  interesting  example  of  the  utility  of  beam  quality  measure¬ 
ments,  we  have  very  recently  begun  a  series  of  measurements  of  the  laser  beam  qual¬ 
ity  and  beam  propagation  factors  in  both  transverse  directions  simultaneously  for  a 
variety  of  semiconductor  diode  lasers,  including  both  narrow-stripe  single-mode  and 
wide-stripe  multimode  diode  lasers.  These  metisurements  are  being  made  using  both 
the  Coherent  beam  quality  measuring  instrument  and  our  own  CCD-camera-based 
beam  quality  meter.  Preliminary  indications  are  that  measurements  of  this  type  can 
bring  out  very  interesting  new  information  concerning  the  modal  properties  and  oscil¬ 
lation  behavior  of  diode  lasers,  including  correlations  between  the  “1  'nk”  phenomena 
that  are  often  observed  in  such  diode  lasers  and  discontinuities  in  the  beam  quality 
parameters.  In  preliminary  tests  of  this  concept  we  have  also  discovered  that  the  beam 
quality  and  other  beam  properties  of  wide-stripe  diode  lasers  can  be  measured  even  in 
the  sub-threshold  ASE  region,  and  the  variation  in  all  of  the  beam  parameters  followed 
continuously  through  the  threshold  region.  Continued  work  in  this  area  is  presently 
being  carried  out. 
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V.  MEASUREMENT  OF  EXCESS  QUANTUM  NOISE  EFFECTS 

By  way  of  background  for  this  section,  in  1979  K.  Petermann  in  Germany  pub¬ 
lished  a  theoretical  analysis  predicting  the  existence  of  an  "excess  spontaneous  emission 
factor”  or  excess  quantum  noise  emission  effect  in  gain-guided  semiconductor  diode 
lasers.  In  1989  under  earlier  AFOSR  support  we  then  published  several  theoretical 
analyses  showing  that  this  kind  of  excess  emission  would  also  occur  and  would  lead 
to  a  sizable  enhancement  of  the  fundamental  quantum  noise  fluctuations  in  any  kind 
of  laser  oscillator  using  an  unstable  optical  resonator.  Our  analysis,  aided  by  theoret¬ 
ical  calculations  using  the  P.A.RAXLA.  software  described  in  Section  III.  showed  that 
the  excess  noise  factors  in  a  typical  unstable  resonator  laser  design  could  range  from 
several  hundred  to  several  thousand,  according  to  our  calculations. 

These  quantum  fluctuations  in  turn  will  represent  the  ultimate  lower  limit  on  the 
amplitude  and  frequency  stability  of  any  laser  device  using  such  a  resonator.  Environ¬ 
mental  or  technical  noise  is  usually  predominant  over  quantum  noise  effects  in  practical 
laser  devices.  Current  developments  in  laser  technology,  however,  are  leading  to  new 
laser  devices  with  extremely  small  noise  fluctuations  and  spectral  widths,  approaching 
the  quantum  limit;  and  such  ultrastable  laser  detdees  are  now  becoming  of  practi¬ 
ced  importance  in  ultrahigh-performance  coherent  optietd  communications  systems,  in 
laser  radars,  and  potentially  also  in  gravity  wave  experiments.  It  becomes  important, 
therefore,  to  verify  experimentally  the  existence  of  this  excess  noise  emission  and  the 
correctness  of  our  theoretical  predictions,  not  only  for  their  fundamental  interest  as 
an  unusual  consequence  of  quantum  electronics,  but  also  in  order  to  set  the  ultimate 
quantum  limits  on  the  performance  of  ultrastable  laser  devices. 

We  have  therefore  been  carrying  out  as  part  of  our  current  work  a  project  to  make 
a  careful  and  definitive  measurement  of  this  excess  noise  factor  in  an  unstable-resonator 
laser.  During  the  previous  and  current  reporting  periods  we  designed  a  miniaturized 
monolithic  diode-pumped  Nd.YAG  laser  system,  employing  diode  pumping  for  pump 
stability;  a  monolithic  rod  structure  for  mechanical  and  acoustic  stability;  and  nitrogen 
cooling  of  the  miniature  YAG  rod  to  narrow  the  laser  linewidth  and  thus  obtain  the 
increased  gain  needed  for  the  unstable  cavity.  During  the  current  reporting  period  we 
have  now  built  and  operated  this  laser  system,  using  a  100  mW  single-mode  diode 
Itiser  pump  and  a  3-mm  long  Nd;Y.4G  rod  which  is  enclosed  inside  a  liquid  nitrogen 
cryogenic  refrigerator.  The  unstable  resonato’  is  ground  directly  onto  the  monolithic 
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rod,  and  the  diameter  of  the  unstable  mode  is  controlled  in  this  initial  experiment 
by  the  inverted  region  created  by  the  focused  diode  laser  pump  beam.  Much  of  the 
experimental  work  up  to  this  {)oint  is  summarized  in  the  conference  proceedings; 

P.  L  Mussche  and  A.  E.  Siegman.  “Enhanced  Schawlow- Townes 
linewidth  in  lasers  with  nonorthogonal  transverse  eigenmodes"  in 
Laser  Noise,  R^jarshi  Roy,  ed..  Proc.  SPIE  1376  (November  1990). 
[Proceedings  of  SPIE  Symposium  on  Laser  Nose,  OPTCON  90, 
November  1990.] 

A  reprint  of  this  report  is  attached  as  Appendix  I.  Since  then  we  have  also  designed 
and  assembled  a  Pound-Drever  cavity  stabilization  system  which  is  used  to  measure 
the  quantum  frequency  fluctuations  of  our  experimental  laser.  In  this  system  a  slow 
feedback  loop  is  used  to  lock  the  reference  cavity  to  the  laser  frequency,  so  that  the 
Pound-Drever  system  tracks  slow  drifts  in  the  laser  frequency.  At  the  same  time,  faster 
fluctuations  which  are  above  the  cutoff  frequency  of  the  feedback  system  are  measured 
in  the  Poimd-Drever  output  signal  and  used  to  determine  the  quantum  noise  spectrum 
of  the  laser  output.  This  system  is  now  in  operation  and  working  well.  The  shot  noise 
in  the  measurement  system  has  been  established  using  a  stable-cavity  diode-pumped 
YAG  laser  with  a  much  lower  Schawlow-Townes  noise  level,  and  the  enhanced  noise  in 
our  unstable  laser  then  measured  with  respect  to  this. 

With  the  present  system  we  have  at  this  point  been  able  to  obtain  a  tentative 
confirmation  of  the  expected  excess  aoi.=^e  level  in  our  unstable  laser;  and  the  excess 
noise  results  we  have  obtained  to  date  on  this  configuration  appear  to  be  in  reasonable 
agreement  with  theory.  We  are  attempting  now  to  fabricate  several  similar  monolithic 
unstable- resonator  YAG  Icisers  with  conventional  hard-edged  mirrors  which  will  rep¬ 
resent  a  wider  range  of  unstable-resonator  parameters,  in  order  to  make  careful  and 
definitive  measurements  of  the  excess  noise  factors  in  these  rods  and  thereby  confirm 
in  a  definitive  fashion  the  excess  noise  figures  we  have  pr  ^dieted.  These  results  should 
then  provide  a  firm  foundation  for  future  design  and  evaluation  of  ultrastable  or  ul¬ 
traquiet  laser  oscillator  and  laser  system  designs,  as  well  as  confirming  the  unusual 
<iuantum  aspects  of  the  predicted  e.xcess  spontaneous  emission. 


VI.  ULTRASHORT  PULSES  IN  THE  INFRARED 


Finally,  in  previous  work  under  this  contract  we  have  summarized  the  poten¬ 
tial  advantages  and  capabilities  of  synch-pumped  and  mode-locked  optical  parametric 
oscillators  (OPOs)  pumped  by  mode-locked  NdrYAG  lasers.  Mode-locked  OPOs  pro¬ 
vide  a  promising  method  for  generating  very  short  pulses  in  the  infrared;  and  if  these 
OPOs  can  be  pumped  by  diode-pumped  Nd:YAG  lasers,  one  has  the  potential  for  a 
comparatively  efficient  all-solid-state  sj'stem  for  IR  pulse  generation. 

In  our  work  during  the  current  period  we  have  constructed  and  operated  two 
such  synch-pumped  and  mode-locked  OPO  systems,  making  use  in  the  first  case  of  a 
LiNbOa  crystal  pumped  at  the  doubled- YAG  wavelength  of  532  nm  and  oscillating  at 
signal  and  idler  wavelengths  (jf  ~SS0  nm  and  ~1350  nm,  and  in  the  second  case  of  a 
KTP  crystal  pumped  at  the  fundamental  1064  nm  YAG  wavelength  with  signal  and 
idler  at  ~1.5  microns  and  ~3.5  microns.  Both  of  these  systems  can  now  be  operated 
reliably  well  above  threshold  with  the  ^'AG  pump  laser  operating  in  a  mode-locked 
plus  repetitively  Q-switched  mode  (at  a  1  kHz  repetition  rate);  and  both  are  clearly 
mode-locking  and  generating  short  IR  pulses.  We  are  now  proceeding  with  efforts  to 
measure  the  exact  pulsewidths  and  other  detailed  behavior  of  these  systems. 

Our  future  objectives  on  this  project  will  include  making  better  and  more  detailed 
metisurements  of  the  operation  of  these  mode-locked  OPOs  in  the  repetitive  Q-switched 
mode,  including  determination  of  pulse  width  at  the  IR  frequencies,  especially  for 
the  KTP  oscillator  in  the  middle  IR.  In  addition,  we  will  attempt  to  optimize  the 
performance  and  operation  of  these  OPOs,  including  reducing  the  pump  threshold 
and  increasing  our  Nd:YAG  laser  performance  so  that  one  or  both  of  these  OPOs — 
with  the  KTP  version  again  taking  precedence — can  be  operated  as  a  reliable  CW 
mode-locked  system. 

We  are  also  considering  an  effort  to  take  advantage  of  a  combination  of  self  phase 
modulation  together  with  the  ncwdy  emerging  concept  of  "Kerr  lens  mode  locking", 
to  obtain  significant  pulse-shortening  in  our  mode-locked  Nd:Y.4G  pump  laser,  both 
for  its  utility  in  our  mode-locked  OPO  experiments,  and  as  a  potentially  simple  and 
inexpensive  way  of  retrofitting  a  large  number  of  existing  Y.4G  lasers  to  obtain  sub¬ 
stantially  shorter  mode-locketl  i)ulses. 
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VII.  HONORS  AND  AWARDS 


During  the  current  reporting  period  Professor  Siegman  presented  the  R.V.  Pole 
Memorial  Lecture,  “Laser  Resonators;  Retrospective  and  Prospective’’ .  at  the  1990 
CLEO/QELS  Meeting  in  Anaheim,  CA,  May  21,  1990. 

Professor  Siegman  also  presented  a  George  Eastman  Invited  Lecture  on  “Defining 
and  Measuring  Laser  Beam  Quality,  "  to  the  Optical  Society  of  America,  Washington 
DC,  April  18,  1990. 

In  November  of  this  year  Professor  Siegman  will  receive  the  1991  Schawlow  Award 
of  the  Laser  Institute  of  America  and  will  be  the  Honored  Speaker  at  the  LIA’s  Annual 
ICALEO  Meeting  in  San  Jose. 
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AbsIracl—V/e  present  a  new  formalism  to  describe  beam  propaga¬ 
tion  in  paraxial  optical  systems  with  dispersive  elements,  including  both 
spatial  and  temporal  variations  in  the  propagating  signal.  This  new 
formalism  makes  use  of  4  x  4  “ray-pulse”  matrices  which  take  ac¬ 
count  of  dispersive  effects  up  to  quadratic  phases  in  both  spatial  co¬ 
ordinates  (as  in  the  usual  paraxial  ABCD  matrix  approach)  and  in  the 
temporal  domain.  We  show  how  to  use  these  matrices  to  write  a  space- 
time  integral  analogous  to  a  generalized  Huygens  integral,  and  derive 
propagation  laws  for  Gaussian  ray  pulses  which  are  space-  and  time- 
varying  analogs  of  the  conventional  results  for  Gaussian  beams.  This 
new  formalism  should  be  very  useful  in  analyzing  dispersive  optical 
systems  such  as  prism  beam  expanders,  femtosecond  pulse  compres¬ 
sion  systems,  and  dispersive  mode-locked  laser  cavities. 

I.  Introduction 

Martinez  [l],  [2]  has  recently  discussed  the  use  of 
3x3  matrices  to  model  polychromatic  beam  prop¬ 
agation  through  paraxial  optical  systems  with  dispersive 
elements.  His  treatment  starts  by  defining  an  optical  axis 
as  the  path  a  midband  reference  ray  takes  through  the  sys¬ 
tem.  He  then  constructs  3x3  matrices  of  the  same  form 
given  in  Siegman  [3)  which  model  the  path  of  slightly 
displaced  and  slightly  detuned  rays.  In  particular,  the  E 
and  E  elements  of  the  matrix  are  frequency  dependent, 
allowing  different  wavelengths  to  traverse  different  paths 
through  the  system. 

Several  aspects  of  this  treatment  can  be  improved  by 
expanding  the  formalism.  First,  the  Martinez  formulation 
has  ray  matrix  elements  which  depend  upon  the  frequency 
of  the  incident  wave,  whereas  one  would  like  the  vector 
that  describes  the  ray  to  carry  all  of  the  information  about 
the  light  and  the  matrix  that  represents  the  optical  system 
to  have  no  dependence  on  the  ray.  Second,  the  properties 
of  the  optical  system  are  not  fully  described  by  the  3  x  3 
matrix,  since  several  additional  phases  must  be  carried 
along  to  get  a  complete  description  of  the  optical  system. 

In  this  paper,  we  show  how  a  formalism  making  use  of 
4x4  matrices  can  be  used  to  model  time-  and  space- 
varying  wave  propagation  in  dispersive  paraxial  systems 
in  a  canonical  way  that  overcomes  these  difficulties.  In 
particular,  the  optical  system  is  completely  modeled  by 

Manuscript  received  July  5.  1989;  revised  January  3,  1990.  This  work 
was  supported  in  pan  by  the  Air  Force  Office  of  Scientific  Research  and 
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its  ray-pulse  matrix,  while  the  vector  that  represents  the 
ray  contains  all  the  information  about  both  its  temporal 
frequency  and  its  position. 

II.  The  Ray-Pulse  Matrix  Formalism 

We  begin  by  constructing  an  optic  axis  which  corre¬ 
sponds  to  the  path  of  a  spatially  transform-limited,  pulsed, 
midband  reference  beam  as  it  goes  through  the  optical 
system.  At  each  point  along  this  reference  beam,  we  erect 
transverse  coordinates  in  the  plane  perpendicular  to  the 
reference  ray.  For  convenience,  we  place  the  origin  of 
these  coordinates  at  the  reference  beam,  as  shown  in  Fig. 

1 .  This  beam  is  also  to  be  thought  of  as  being  a  temporally 
transform-limited  pulse,  so  that  it  can  simultaneously  have 
a  well-defined  center  frequency  (the  reference  frequency 
or  central  frequency  for  the  system)  and  can  mark  a  well- 
defined  time  of  arrival  at  each  transverse  plane.  For  con¬ 
venience,  we  set  the  clocks  in  each  transverse  plane  to 
zero  at  the  time  of  arrival  of  the  reference  pulse. 

We  can  now  imagine  firing  another  transform-limited 
pulse  at  a  different  starting  time  and  with  a  different  car¬ 
rier  frequency  into  the  system  at  a  slightly  different  po¬ 
sition  and  slope.  We  will  denote  the  changes  in  these  four 
quantities  as  and  0|n,  respectively.  In  order  to 

simplify  later  results,  we  assume  that  each  reference  plane 
is  surrounded  by  a  slab  of  infinitesimal  thickness  of  ma¬ 
terial  of  unit  index  of  refraction.  Equivalently,  the  0’s 
should  be  viewed  as  the  slope  normalized  by  the  index  of 
refraction.  This  normalization  avoids  the  need  to  have 
matrices  that  describe  flat  dielectric  interfaces.  Also,  we 
use  hertzian  frequency  /  rather  than  radian  frequency  w  as 
this  reduces  the  number  of  tt's  that  need  to  be  carried 
along. 

As  this  new  pulse  crosses  each  of  the  transverse  planes, 
it  defines  a  new  arrival  time,  frequency,  position,  and 
slope  relative  to  the  reference  pulse.  (The  time  at  each 
plane  is  really  a  time  offset  from  the  reference  time  at  that 
plane.)  We  know  that  the  coordinates  at  the  output  plane 
must  be  a  function  of  the  input  parameters,  i.e.,  that 
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Fig.  1.  Optical  system  with  reference  planes  (R.P.  1-R.P.7)  and  clocks 
constructed  to  measure  path  and  time  of  flight  of  arbitrary  pulse.  Note 
that  all  clocks  in  R.P.  1-R.P.5  show  positive  times,  but  those  in  R.P. 6 
and  R.P. 7  are  negative  and  will  reach  zero  when  reference  pulse  arrives. 
Each  R.P.  is  perpendicular  to  direction  of  propagation  of  reference  pulse, 
and  has  its  origin  where  reference  pulse  goes  through  R.P. 


Also,  by  the  definition  of  the  reference  ray,  we  know  that 
SF  { 0 }  =0.  For  systems  of  interest  the  function  fF  is 
highly  differentiable,  so  that  it  makes  sense  to  replace  (1) 
by  the  first  nonzero  term  in  its  power  series.  This  yields 
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Most  of  these  matrix  elements  can  be  put  in  simpler  terms. 
Consider  first  the  upper  left  2x2  block  in  the  matrix. 
These  elements  determine  the  changes  in  output  position 
and  slope  caused  by  input  changes  in  position  and  slope, 
with  input  frequency  and  time  unchanged,  and  are  thus 
just  the  conventional  ABCD  matrix  elements  for  the  sys¬ 
tem.  (This  is  not  entirely  obvious,  as  the  conventional 
ABCD  matrix  is  defined  in  terms  of  the  propagation  of  a 
monochromatic  wave,  while  the  elements  of  the  general¬ 
ized  ray  matrix  are  defined  in  terms  of  the  propagation  of 
a  group.  A  rigorous  proof  of  this  equivalence  will  be  given 
below  at  the  end  of  the  section  on  the  generalized  Huy¬ 
gens  integral.) 

Similarly,  looking  at  the  bottom  row  and  knowing  that 
in  all  cases for  a  time  invariant  linear  system, 
we  find  that  this  row  must  be  0,  0,  0,  1 .  Finally,  the  third 
column,  which  determines  how  the  output  pulse  depends 
on  the  input  time  must  be  0,  0,  1,0  due  to  the  time  in¬ 
variance  of  the  system.  Putting  all  of  these  simplifications 


in  place,  we  have 


A 

B 

0 

3j:out 

3fin 

X 

33„ui 

X 

C 

D 

0 

3 

3f|n 

9 

= 

X 

t 

^^oul 

3tout 

1 

3tout 

t 

J  _ 

out 

^Xjf, 

33.„ 

3f,„ 

f 

0 

0 

0 

1 

(3) 


For  convenience,  we  name  the  remaining  derivatives  E 
thru  /,  so  that  a  generalized  4  x  4  ray  matrix  or  ray-pulse 
matrix  has  the  form 
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Each  of  the  symbolic  entries  in  this  matrix  can,  in  fact, 
take  on  any  value,  although  the  nine  entries  are  not  all 
independent.  One  of  the  relations  among  them  is  the  fa¬ 
miliar  AD  -  BC  =  1 ,  while  two  other  relations  involve 
the  dispersive  elements  and  have  no  analog  in  ABCD  ma¬ 
trices.  The  relations  will  be  discussed  in-depth  in  the  sec¬ 
tion  on  the  analog  of  the  Huygens  integral,  where  the  form 
of  the  extra  relations  is  easy  to  discover;  however,  it  is 
easy  to  see  why  AD  -  BC  =  1  without  the  use  of  the 
integral  formalism.  If  we  imagine  placing  blackbodies  at 
the  same  temperature  on  both  ends  of  the  optical  system, 
each  emits  a  certain  amount  of  energy  per  unit  time  per 
unit  frequency  per  unit  area  into  a  unit  angle.  Because  the 
matrix  in  (4)  gives  the  output  time,  frequency,  position, 
and  slope  as  functions  of  the  input,  its  determinant  is  the 
magnification  of  the  phase  space  volume  that  the  radiation 
is  in.  But  by  the  second  law  of  thermodynamics,  there  can 
be  no  net  power  transferred,  and  thus  no  change  in  radia¬ 
tion  density,  so  that  the  determinant  must  be  1.  Expand¬ 
ing  out  the  determinant  then  yields  AD  -  BC  =  1 . 


III.  Ray-Pulse  Matrices  for  some  Common 
Elements 

As  an  example  of  this  formalism,  let  us  compute  the 
ray-pulse  matrix  for  a  grating.  The  reference  planes  are 
taken  perpendicular  to  the  incident  and  outgoing  parts  of 
the  reference  ray,  and  are  chosen  for  convenience  to  meet 
at  the  point  of  reflection,  as  shown  in  Fig.  2(a).  Construc- 
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Fig.  2.  Ray-pulse  matrix  for  grating  is  computed  from  difference  in  be¬ 
havior  of  reference  ray  (a),  spatially  offset  ray  (b),  tilted  ray  (c).  and 
detuned  ray  (d). 


tion  of  the  ray-pulse  matrix  begins  by  considering  an  in¬ 
cident  wave  parallel  to  the  reference  ray,  but  separated 
from  it  by  x  units  as  shown  in  Fig.  2(b).  This  new  ray 
strikes  the  grating  jr/sin  units  away  from  the  reference, 
and  the  reflected  ray  can  be  extended  backwards  to  the 
output  plane  where  its  height  is  -jc  sin  </)/sin  Thus,  A 
=  -sin  </)/sin  tt.  From  the  grating  law,  it  is  clear  that  the 
output  is  parallel  to  the  reference  from  which  we  conclude 
that  C  =  0.  Also,  the  extra  distance  traveled  by  the  offset 
ray  is  just  PQ  -  QR  so  that  the  pulse  will  arrive  .v(cos  ^ 
-  cos  <f>)/(c  sin  1^)  later,  so  that  G  =  (cos  i/-  -  cos  </>)/ 
(c  sin  ^) . 

We  now  consider  a  ray  that  meets  the  input  plane  at  the 
same  place  as  the  reference  ray,  but  that  has  a  different 
slope  as  shown  in  Fig.  2(c).  It  is  clear  that  the  output 
position  has  not  moved,  and  thus  that  B  =  0.  From  dif¬ 
ferentiating  the  grating  equation,  (/o  +  /)  (cos  <t>  -  cos 
fir)  =  c/Xgratmg  where  /o  is  the  frequency  of  the  reference 
ray,  we  can  deduce  that  -sin  \lr6\lr  +  sin  <f)5<f)  =  0  so  that 
D  =  -sin  ^/sin  <f>  because  the  input  angle  is  -5f{r.  We 
could  also  have  obtained  this  element  from  the  previously 
obtained  values  for  A  and  B  plus  the  relationship  AD  — 
BC  =  1 .  Also,  the  delay  of  the  tilted  ray  is  seen  to  be 
zero  so  that  H  =  0. 

Finally,  consider  a  pulse  that  lies  on  top  of  the  refer¬ 
ence  pulse,  but  has  a  different  center  frequency,  as  in  Fig. 
2(d).  Just  as  in  the  case  of  the  tilted  input  wave,  it  is  clear 


TABLE  1 

Generalized  Ray  Matrice.s  for  Various  Optical  Components 
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that  the  output  time  and  position  are  identical  with  the 
reference,  and  thus  that  £  =  /  =  0.  To  obtain  the  last 
element,  differentiate  the  grating  equation  with  respect  to 
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'"to  obtain / ( cos  ^  -  cos  <t>)  —  fo  sin  060  =  0  from  which 
ive  conclude  that  F  =  (cos  0  —  cos  0)/(/o  sin  0)  where 
^0  is  the  frequency  of  the  reference  ray . 

Using  the  method  of  the  preceding  paragraphs,  explicit 
lvalues  of  the  ray-pulse  matrices  for  other  common  dis¬ 
persive  systems  have  also  been  calculated,  as  summarized 
in  Table  1.  It  is  interesting  to  note  that  the  matrices  for  a 
grating  and  for  a  prism  used  at  its  tip  are  formally  iden¬ 
tical  and,  in  fact,  the  only  difference  (to  the  order  of  ex¬ 
pansion  treated  here)  is  that  the  grating  has  much  larger 
dispersive  terms.  Because  of  this,  we  will  explicitly  treat 
optical  systems  that  are  made  with  prisms,  but  the  reader 
should  remember  that  any  of  these  prisms  can  be  replaced 
with  gratings.  We  consider  now  several  uses  of  these  4 
by  4  ray-pulse  matrices. 

IV.  Multiple-Prism  Systems 

These  ray-pulse  matrices  can  be  used  to  simplify  the 
analysis  of  a  variety  of  useful  dispersive  propagation  sys¬ 
tems.  As  a  first  example  of  a  more  complicated  but  useful 
system,  let  us  consider  a  prism  assembly  similar  to  that 
used  in  femtosecond  laser  cavities,  as  discussed  by  Fork 
[4].  We  assume  that  the  reference  beam  passes  through 
the  apexes  of  all  four  prisms;  that  the  input  face  of  the 
first  prism  is  parallel  to  the  output  face  of  the  second 
prism;  that  all  of  the  prism  apex  angles  are  the  same;  and 
that  the  second  pair  of  prisms  is  a  mirror  image  of  the  first 
pair,  as  indicated  in  Fig.  3(a).  These  specific  conditions 
ensure  that  the  rays  at  the  input,  output,  and  between  the 
second  and  third  prisms  are  parallel  independent  of  the 
index  of  refraction  and  dispersion  of  the  glass  in  the 
prisms.  Under  these  conditions,  the  first  prism  has  a  gen¬ 
eralized  ray  matrix  given  by 


Mirror 

Symmetry 

Plane 


M4.  D4  ^  MS,  D5 


Ml.  01 


M2.  02  M3. 03  \  " 

M6.  06  , 

\M7.  07 

(C)  ^ 

Fig.  3.  Various  systems  that  can  be  modeled  with  generalized  ray  matri¬ 
ces.  (a)  Multiple-prism  system  used  to  obtain  positive  group  velocity 
dispersion  in  femtosecond  laser  cavities,  (b)  Variation  on  frequently  used 
for  short  pulse  amplification,  (c)  General  prism  beam  expander  whose 
ith  prism  has  magnification  M, . 


ond  and  third  is  L,  the  overall  system  matrix  is  found  to 
be 


\/M  0  0  0 

0  Af  0  -MD 

-D/)xo  0  1  0 

0  0  0  1 

Explicit  values  of  M  and  D  in  terms  of  the  physical  pa¬ 
rameters  of  the  prisms  are  of  little  importance  here,  but 
can  be  found  in  Table  I.  The  matrices  of  the  third  and 
fourth  prisms  are  respectively  identical  to  (5)  and  (6)  but 
with  D  replaced  by  —D.  Assuming  that  the  space  between 
the  first  and  second  prisms  is  L.  amd  that  between  the  sec- 


This  result  says  that  rays  of  different  wavelength  that  come 
into  the  sy.stem  at  the  same  position  and  slope  will  go  out 
of  the  system  at  the  same  position  and  slope,  but  at  dif¬ 
ferent  times.  Indeed,  this  matrix  is  identical  to  that  for  a 
slab  of  material  with  a  positive  group  velocity  dispersion, 
independent  of  the  sign  of  D,  so  that  the  cascaded  prism 
assembly  provides  a  method  of  achieving  positive  group 
velocity  dispersion  with  no  angular  dispersion. 

This  approach  has  also  been  applied  to  the  specific  case 
discussed  by  Fork,  in  which  the  reference  ray  does  not 
pass  through  the  apexes  of  the  second  and  third  prisms, 
and  the  prisms  have  been  cut  so  that  they  can  be  simul¬ 
taneously  used  at  Brewster's  angle  and  Af  =  1 .  The  result 
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obtained  in  this  way  forc/r/J/are  in  exact  agreement  with 
that  found  by  Fork  by  a  geometric  construction. 

As  a  simple  variation,  the  second  pair  of  prisms  can  be 
arranged  as  in  Fig.  3(b),  and  a  telescope  of  magnification 
- 1  can  be  used  to  fill  the  space  between  the  third  and 
fourth  prisms.  This  yields  an  overall  system  matrix  of 

-1  LOO 
0-100 
0  0  10 

0  0  0  1 

which  results  in  the  beam  being  turned  upside  down  and 
propagated  backwards  by  a  distance  L.  Assemblies  simi¬ 
lar  to  this  are  in  common  use  in  ultrashort  pulse  amplifi¬ 
cation  schemes.  These  systems  are  gratings  instead  of 
prisms,  and  are  designed  to  temporally  disperse  the  pulse 
before  it  passes  through  the  gain  medium  (located  be¬ 
tween  the  second  and  third  grating)  and  then  reassemble 
the  frequency  components  in  the  last  two  gratings. 

As  a  final  example,  consider  the  design  of  an  achro¬ 
matic  prism  beam  expander  as  shown  in  Fig.  3(c),  where 
one  wants  to  have  the  direction  of  the  output  beam  not 
vary  with  the  frequency  of  the  input  (i.e.,  F  =  0).  Since 
we  are  only  interested  in  this  one  matrix  element,  it  is 
clear  that  we  can  omit  the  matrices  that  correspond  to  the 
free  space  between  the  prisms,  and  that  we  may  also  as¬ 
sume  that  the  beams  go  through  the  tips  of  the  prisms.  A 
straightforward  induction  shows  that  the  overall  matrix 
for  k  pri.sms  is  ju.vi 

W*  0  0  0 

0  1 /M*  0  D* 

0  1  0 

0  0  0  1 

Mt  0  0  0 

0  1//W|  0  D| 

0  1  0 

_0  0  0  1 

n  A/, 

/  - 1 

0 

1  I 

1/X„  X  D,  11 

I  [  -It  ^  I 


where  Mj  is  the  magnification  due  to  the  jth  prism,  and 
is  the  rate  of  change  of  output  angle  with  input  frequency 
in  the  jth  prism.  The  overall  angular  dispersion  is  thus 
=  I  =  /  +  I  M„,  in  agreement  with  the  result  implied 

in  Duarte  and  Piper  [5]  and  given  explicitly  in  Trebino 
|6|.  This  result  implies  that  under  common  conditions  the 
optimal  design  for  an  achromatic  prism  beam  expander 
should  have  all  but  the  last  prism  pointing  in  the  same 
direction,  as  was  pointed  out  by  Trebino. 


V.  Generalized  Huygens  Integral 

It  is  well  known  that  optical  systems  characterized  by 
an  ordinary  ABCD  matrix  are  equivalently  represented  by 
the  monochromatic  Huygens  integral  operator  13] 


!  A 


where  we  have  omitted  a  constant  phase  factor.  This  can 
be  thought  of  as  saying  that  the  optical  path  length 
(eikonal)  through  an  arbitrary  ABCD  system  from  .v^  to 
is  ( 1/2B)  (Ax;„  -  2.r,n.r„u,  -l-  D.r,;„,)  plus  a  constant 
which  we  will  ignore. 

In  this  section,  we  will  derive  an  analog  to  this  result 
for  the  4x4  ray-pulse  matrices.  We  begin  by  specifying 
both  the  position  and  time  that  a  pulse  is  to  arrive  at  the 
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input  and  output  planes.  (This  ability  to  use  two-point 
boundary  conditions  is  a  result  of  the  even  number  of  co¬ 
ordinates,  since  a  system  with  an  odd  number  of  coordi¬ 
nates  will  either  be  underdetermined  or  overdetermined 
when  an  even  number  of  conditions  are  applied.)  The 
specification  of  these  four  coordinates  allows  unique  re¬ 
construction  of  the  other  four  coordinates,  except  in 
pathological  cases  such  as  imaging  between  the  input  and 
output  planes. 

The  integral  kernel  is  now  expected  to  be  of  the  form 
exp  [  -y27rZ,(.r,„,  r„ut)/^o]  where  the  function  L 

is  real.  Expanding  L  in  a  power  series  in  all  four  coordi¬ 
nates  yields  first  a  constant  length  Z-o  that  is  analogous  to 
the  overall  phase  factor  in  the  Huygens  integral  formula¬ 
tion,  which  we  will  ignore  just  as  we  have  done  in  (10). 
Next,  there  are  linear  terms  that  are  zero  because  of  the 
way  we  have  constructed  the  reference  planes.  Again,  this 
is  analogous  to  the  Huygens  formulation  where  there  are 
no  linear  terms  in  the  exponential  when  treating  systems 
in  which  ail  elements  are  on  axis  but,  when  elements  are 
off  axis,  one  must  use  ABCDEF  matrices,  and  the  kernel 
does  indeed  contain  linear  terms  as  given  below  in  (20). 
The  next  term  is  a  quadratic  form  in  x,„,  jr„u,,  and 
that  can  be  reduced  to  a  quadratic  form  in  jcip,  .v„u,,  and 
-  r„u,  by  the  time  invariance  of  the  system.  As  in  the 
Huygens  formulation,  we  stop  with  the  quadratic  terms, 
although  there  are  in  fact  higher  order  terms  that  are  of 
importance  when  rays  are  far  off  axis  or  are  wildly  de¬ 
tuned  from  the  reference  ray.  Thus,  we  have  the  follow¬ 
ing  ansatz  for  the  space  and  time-varying  analog  of  the 
Huygens  integral: 


r  r 

1 

t' 

a 

V  j  J  dx,„  dt,„  U(.r,„,  f,J  exp 

/«  &  7\  A.n  \  ' 

■1/3  b  Ell  .r„„,  I  (11) 

n  f  S'/  \f,n  -  U|/  - 

where  X(,  is  the  wavelength  of  the  reference  ray. 


A  simple  way  to  establish  the  correspondence  between 
the  elements  of  the  generalized  ray  matrix  and  the  un¬ 
known  parameters  a,  (3,  •  •  ■  ,  r/  is  to  take  an  incident 
wave  which  has  its  center  at  position  .r,  time  /,  mean  slope 
6,  and  mean  frequency  f  and  perform  the  integration  in¬ 
dicated  in  (1 1).  Using  the  input  wave 


exp 


2ir 


-  xY  (r,„  - /)>- 


K 


-  X)  +  jf(t,„  -  0 


(12) 


yields  an  output  wave  which  is  a  constant  involving  rj 
times  the  exponential  of  a  second-degree  polynomial  in 
jr„u,  and  with  coefficients  that  depend  upon  a,  (3, 

•  •  •  ,  S'  and  the  coordinates  of  the  input  wave.  In  order 
to  find  the  position,  slope,  frequency,  and  time  of  the  out¬ 
put  pulse,  we  must  separate  the  wave  into  the  product  of 
an  envelope  and  a  phase.  This  is  easily  done  as  the  poly¬ 
nomial  in  the  exponential  is  easily  split  into  real  and 
imaginary  parts,  the  exponential  of  the  real  part  yielding 
the  envelope,  and  the  exponential  of  the  imaginary  part 
being  the  phase.  The  position  and  time  of  the  output  pulse 
can  then  be  located  by  setting  the  derivatives  with  respect 
to  3r„u,  and  fom  of  the  envelope  to  zero  and  solving  forT:„u, 
and  r„u,. 

Alternately,  we  could  evaluate  jj  licoui 
-Vou.  I  I  VJJ  dx,,^,  1  I  -  as  the  mean  position  of 

the  output  pulse,  and  use  JJ  dx„^,  f„„,l  V!I  dx,,^, 

di„ut  I  I  “  as  the  mean  time  for  the  output  pulse.  How¬ 
ever,  a  basic  theorem  for  Gaussian  distributions  tells  us 
that  these  are  identical  to  the  position  and  time  of  the 
maximum  amplitude  of  the  pulse.  Similarly,  the  fre¬ 
quency  and  slope  of  the  output  ray  may  be  found  by  te¬ 
dious  integration,  or  by  evaluating  the  derivative  of  the 
imaginary  part  of  the  polynomial  at  the  position  and  time 
found  Just  above.  The  result  of  this  can  be  summarized  as 
the  generalized  ray  matrix 
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Intermediate  steps  in  the  derivation  of  this  result  have  been 
suppressed  as  they  are  exceedingly  complicated,  and  re¬ 
quired  the  use  of  MACSYMA.  This  result  can  be  inverted 
to  give  a,  P,  ■  •  •  ,  f  as  functions  of  A,  B,  ■  •  ■  ,  /.  The 
results  are 

\oA(EH  -  BI)  -I-  E- 
\iB{EH  -  BI) 


may  be  taken  to  be 

AD  -  BC  =  I 
BE  -  ED  =  \)H 

AF-EC  =  \)G  (15) 

or  they  may  be  taken  to  be 

/ID  -  BC  =  1 


EH  -  BI 
E 

EH  -  BI 

D(EH  -  BI)  +  Xo//- 
B(EH  -  BI) 

EH  -  BI 
X„B 

EH  -  BI 


(14) 


Finally,  rj  can  be  determined  up  to  a  phase  by  insisting 
that  energy  is  conserved.  Doing  so  yields  ri  = 

'J\^(EH  -  BI)  where  the  phase  of  rt  has  been  chosen 
to  make  the  output  phase  identical  to  the  input  in  the  limit 
of  an  identity  ray-pulse  matrix.  With  these  identifications, 
(11)  becomes  a  generalized  time  and  space-varying  Huy¬ 
gens  kernel,  analogous  to  the  space-varying  integral  of 
(10),  and  dependent  only  on  the  A  through  /  matrix  ele¬ 
ments. 

In  the  expressions  for  a,  /3,  •  ■  •  ,  tj  we  find  that  the 
expression  EH  —  BI  occurs  repeatedly  in  the  denomina¬ 
tor.  Several  of  the  systems  that  we  have  already  discussed 
make  this  quantity  zero,  so  we  will  describe  how  to  han¬ 
dle  this  case.  There  are  two  distinct  possibilities;  if  at  leas! 
one  of  E,  H,  B,  or  I  is  nonzero,  the  kernel  collapses  into 
a  one-dimensional  delta  function,  which  cancels  one  in¬ 
tegral  in  (1 1)  and  the  other  integration  remains  nontrivial. 
If,  on  the  other  hand,  all  of  £,  H  S,  and  I  are  zero,  we 
obtain  two  delta  functions  that  renders  both  integrations 
trivial.  Both  of  these  cases  are  handled  by  making  small 
perturbations  in  the  matrix  elements,  and  then  letting  these 
changes  go  to  zero.  For  the  case  EH  -  BI  -  0,  but  B  ^ 
0,  we  obtain  the  kernel  V/ /XqB  exp  [(  -jir/\)B)  (/tcjn 

2xjnJc„yt  4“  Dx,,y()]  6(rin  Cui  +  4“  Cxj^/XoB). 

Formulas  of  a  similar  type  exist  for  E  ^  0,  H  ^  0,  and 
1^0.  For  the  case  where  B  =  £  =  //  =  /  =  0the 
propagator  is  just  exp  (  -jitCx^^,/\QA]  6(x,n  - 

)  ^(^in  Cut  )  . 

Since  the  propagator  has  only  six  parameters  that  influ¬ 
ence  the  output  coordinates  and  the  generalized  ray  matrix 
has  nine  variable  elements,  we  know  that  there  are  three 
relations  between  these  matrix  elements.  These  relations 


Xo(BG  -  AH)  -  E 

Xo(DG  -  CD)  =  F  (16) 

depending  upon  whether  one  wants  to  think  of  H  and  G 
or  E  and  F  as  being  more  fundamental.  These  constraints 
may  also  be  written  in  the  matrix  form 
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(17) 


Quick  inspection  reveals  that  all  of  the  matrices  in  Table 
I  obey  these  relations.  Also,  it  is  easy  to  show  that  if  a 
series  of  ray-pulse  matrices  satisfy  (15),  (16),  or  (17), 
then  .so  does  their  product,  from  which  we  conclude  that 
(15).  (1 6),  and  (1 7)  are  in  fact  satisfied  for  any  system 
made  up  of  the  components  in  Table  1. 

It  was  mentioned  earlier  that  the  A,  B,  C,  and  D  ele¬ 
ments  of  a  generalized  ray-pulse  matrix  are  identical  to 
the  elements  of  a  conventional  ray  matrix,  despite  the  fact 
that  the  ray-pulse  matrix  is  defined  in  terms  of  a  temporal 
group  while  the  conventional  ray  matrix  is  defined  in  terms 
of  a  monochromatic  wave.  There  are  two  distinct  argu¬ 
ments  that  can  be  made  to  establish  this  fact.  First,  an 
intuitive  argument:  note  that  a  temporal  group  can  be  bro¬ 
ken  down  into  monochromatic  waves,  and  that  each  of 
these  can  be  propagated  through  the  system  by  the  use  of 
a  conventional  /4BCD£F  matrix.  The  matrix  used  for  each 
frequency  component  will  be  slightly  diflferent  if  the  sys¬ 
tem  contains  any  dispersive  elements,  which  means  that 
the  various  monochromatic  waves  will  emerge  at  loca¬ 
tions  and  slopes  centered  around  that  of  the  monochro¬ 
matic  wave  with  the  central  frequency.  Finally,  noting 
that  the  energy  per  unit  area  per  unit  angle  of  the  mono¬ 
chromatic  waves  is  additive,  we  find  that  the  temporal 
group  emerges  at  the  location  and  slope  of  its  central  fre- 
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uency  component  that  implies  the  equivalence  of  the  A, 
(,  C,  and  D  elements  of  the  two  types  of  matrix. 

A  second  and  more  rigorous  argument  is  based  upon  the 
luygens  integral  formulation.  We  begin  by  taking  an  ar- 
itrary  generalized  ray  matrix  and  computing  its  gener- 
lized  Huygens  integral  equivalent  as  given  in  (11)  and 
14).  We  then  apply  (11)  to  a  wavefront  C/(Xi„)  with  no 
ime  dependence,  and  perform  the  indicated  time  integral, 
'he  result  is 


^in  ^(-^in) 


exp 


\)B  ^in-*oul  ^^oat) 


(18) 


vhich  is  identical  to  (10)  except  that  the  A,  B,  and  D  ap- 
learing  now  are  elements  of  the  generalized  ray  matrix. 
Thus,  the  generalized  A,  B,  and  D  elements  of  the  ray- 
mlse  matrix  are  identical  to  those  of  the  monochromatic 
ay  matrix.  Finally,  the  C  elements  of  the  two  types 
»f  matrices  are  identical  because  both  satisfy  AD  - 
SC  =  1. 


VI.  Monochromatic  Wave  Propagation 
Now  that  we  have  the  analog  of  Huygens  integral  for 
ay-pulse  matrices,  we  may  obtain  results  for  various  spe- 
:ial  input  wavefronts.  We  begin  with  a  discussion  of  the 
rase  where  the  input  wave  is  monochromatic  at  some  fre- 
luency  /.  In  this  ca.se,  we  will  write  the  incident  wave  as 
[/(.tin)  e^^'^and  perform  the  time  integral.  Doing  so  re¬ 
mits  in  the  output  wave  being  given  by 


exp 


j2^ft  +  {EH  ~  B!)  J  dLr,„  U{xJ 


•  exp 


21- 


vhere  the  eikonal  length  L{x^„.  .tom)  is  given  by 


E{x,„,x^^)  = 


A  -1 

- 1  D2  \x 


+  2 


XoH 


(19a) 


(19b) 


If  the  incident  field  is  known  in  the  form  U{x,f),  the 
Tutput  field  may  be  found  in  the  same  form  by  using  (19) 
for  each  frequency  component.  The  kernel 


—  Vr'" 

^\b[bF-  De) 


(20) 


corresponds  to  an  ABCDEF  matrix  at  wavelength  X  so  that 
the  integral  in  (19)  corresponds  to  the  matrix 


-A  B  fE- 
C  D  fF 
-0  0  1- 


(21) 


plus  terms  that  are  ''quadratic  and  higher  in  the  input 
variables.  These  higher  order  terms  are  not  to  be  trusted 
as  the  entire  theory  can  only  predict  accurately  changes 
in  position  that  are  linear  in  the  input  variables.  This  re¬ 
sult  generalizes  the  argument  made  in  the  last  section 
about  the  correspondence  between  the  A  through  D  ele¬ 
ments  of  the  conventional  and  generalized  ray  matrices, 
in  that  it  also  shows  that  the  E  and  F  elements  of  the  ray- 
pulse  matrix  are  closely  related  to  the  E  and  F  elements 
of  conventional  monochromatic  but  misaligned  ABCDEF 
matrices. 

At  this  point,  the  problem  with  using  only  3x3 
ABCDEF  matrices  to  model  dispersive  systems  also  be¬ 
comes  clear:  the  elements  G,  H,  and  /  are  missing  from 
the  3  X  3  formulation.  The  loss  of  G  and  H  causes  no 
loss  of  information,  as  they  are  computable  from  the  re¬ 
maining  elements  via  (15),  but  /,  which  contains  infor¬ 
mation  about  the  dispersion  on  axis  is  missing  and  is  not 
recoverable.  It  is  important  to  note  that  the  reduction  of  a 
ray-pulse  matrix  to  an  ABCDEF  matrix  for  each  frequency 
ignores  the  phase  factor  in  front  of  the  integral  in  (19), 
and  thus  results  in  the  inability  of  the  reduced  formalism 
to  predict  the  relative  phase  of  various  frequency  com¬ 
ponents. 

A  physical  description  of  this  loss  of  information  may 
be  obtained  from  the  standard  construction  of  the  eikonal. 
For  a  single  frequency,  one  can  develop  the  phase  of 
propagation  through  the  system  by  considering  rays  that 
originate  at  a  point  some  distance  back  from  the  optical 
system,  and  come  to  a  focus  downstream  from  the  sys¬ 
tem.  Since  each  of  these  rays  must  have  the  same  total 
optical  length,  and  the  paths  outside  of  the  system  are 
known,  we  can  subtract  to  get  the  path  in  the  system. 
Thus,  the  system  path  length  involves  the  total  length 
which  may  vary  from  one  frequency  to  the  next.  If  we 
were  to  treat  the  system  with  ABCDEF  matrices,  we  would 
have  to  carry  this  extra  path  or  phase  as  the  auxiliary  func¬ 
tion  Martinez  calls  <t>n-  H  we  make  a  Taylor-series  expan¬ 
sion  of  <>0,  the  constant  term  is  unimportant;  the  linear 
term  is  zero  because  of  the  way  we  have  set  up  the  clocks; 
and  the  quadratic  term  is  given  in  (19).  Thus.  (t>o  is  in¬ 
cluded  in  our  theory,  provided  that  we  neglect  cubic  and 
higher  terms. 


VII.  Generalizf-D  Gaussian  Beam  Propagation 

Another  spjecial  input  wave  that  results  in  interesting 
output  waves  is  a  generalized  Gaussian  pulse.  This  can  be 
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written  as  a  wave  of  the  form 


which  is  a  generalized  form  of  the  usual  Gaussian  beam 
with  beam  parameter  q.  The  minus  sign  in  the  definition 
(22)  is  needed  to  make  the  Q  matrix  have  one  covariant 
and  one  contravariant  index,  as  is  required  by  the  formula 
below.  The  generalized  Gaussian  beam  parameter  C  can 
be  defined  without  the  minus  sign,  but  this  results  in  sev¬ 
eral  appearances  of  the  metric  tensor  below.  Also,  since 
Q  ~ '  defines  a  quadratic  form,  we  require  that  the  xt  and 
lx  elements  of  ^  '  are  negatives  of  each  other. 

Several  other  points  also  need  to  be  made  about  (22). 
First,  the  xx  element  of  Q  does  not  give  the  beam  width 
unless  the  xt  element  is  zero.  Similarly,  the  tt  element 
does  not  determine  the  pulse  width  unless  the  cross  term 
is  zero.  Finally,  note  that  the  intensity  and  phase  fronts 
of  (22)  may  be  obtained  by  taking  the  imaginary  and  real 
parts  of  ^  respectively.  It  is  then  not  only  possible  for 
the  two  real  quadratic  forms  that  result  to  not  have  prin¬ 
cipal  axes  along  .r  and  t.  but  it  is  also  possible  for  them 
to  have  different  principle  axes  that  result  in  pulses  that 
have  phase  fronts  and  intensity  fronts  pointed  in  different 
directions. 

If  the  pulse  (22)  is  propagated  through  the  sy.stem  by 
the  Huygens  integral  (II)  above,  we  find  that 


^IHIt 


£/\.y 

D  F/K\ 
0  I  / 


(23) 


which  is  analogous  to  the  result  q„a,  =  {Aq^^  t  B)/{Cq,„ 
+  D)  that  holds  for  Gaussian  beam:  .  Of  course,  the  an¬ 
alog  is  not  quite  exact,  as  the  expressions  previous  are 
matrices  and  the  order  of  the  numerator  and  denominator 
matters,  while  in  the  conventional  formula  no  such  dis¬ 
tinction  needs  to  be  made.  The  2x2  blocks  in  the  expres¬ 
sion  above  can  be  written  down  by  inspection  by  the  fol¬ 
lowing  means:  divide  the  variable  elements  in  the  last 
column  by  Xo,  and  then  reorder  the  rows  and  columns  of 
the  ray-pulse  matrix  so  that  the  variables  occur  in  the  or¬ 
der  x,  t,  6,  and  /.  The  upper  right  2x2  block  is  the 
coefficient  of  in  the  numerator,  the  upper  left  2x2 
block  is  the  constant  in  the  numerator,  and  similarly  the 
lower  half  of  the  rearranged  matrix  yields  the  denomina¬ 
tor  of  the  expression  for 

We  can  now  try  to  determine  an  eigen-0  for  the  system. 
A  quick  computation  shows  that  there  is  no  input  Q  that 
can  yield  the  same  Q  at  the  output  except  in  degenerate 
cases,  and  then  any  value  will  do  for  The  reason  for 
this  is  that  the  time  part  of  the  problem  behaves  in  a  di¬ 
vergent  fashion:  the  temporal  width  of  any  pulse  expands 


without  bound  as  it  is  repeatedly  passed  through  a  (time 
invariant)  dispersive  system,  although  it  may  narrow  for 
the  first  few  passes.  This  means  that,  except  in  the  case 
of  no  dispersion,  there  cannot  be  any  temporal  eigen- 
pulse,  while  in  the  case  of  zero  dispersion,  any  temporal 
pulse  shape  will  come  through  unaltered.  Returning  to  the 
combined  space  and  time  treatment,  it  is  reasonable  to 
look  for  pulse  hapes  in  which  0,,  and  do  not  change 
upon  propagaiiv,..  'hrough  the  system,  while  Q,,  increases 
by  A  in  each  pass.  Expanding  out  (23)  under  these  as¬ 
sumptions  yields 


(2« 


AQ,,  +  B 
CQ.,  +  D 


Q.  -  GQ,,  -  H/K, 
~  +  D 

CQ„  +  F/\,  _ 


-  GQ„  -  I/\  -f  Q„.  (24) 


Letting  R  =  -JiA  -t-  D)*  -  4,  we  obtain 


A  -  D  ±  R 
2C 


AG  -  F/\^  +  CH  ±R 

~  C{1  +  A  -  D  ±  R) 

AG  -  F/X<,  +  CH  ±  R 
A  =  — //Xf)  —  . 

C{2  +  A  -  D  ±  R) 

|/1G  -  F/Xo  +  CH  ±  R  +  (F/Xo  +  C) 

■{2  +  A-D±R)]  (25) 

as  the  parameters  of  the  quasi-eigenpuise.  Although  the 
result  in  (25)  for  appears  to  be  familiar  from  cavity 
mode  theory,  the  beam  waist  is  not  VAoj5„/’^  unless 
is  zero.  When  the  optical  system  is  designed  to  make 
propagation  separable  into  space  and  time  components. 
(25)  reduces  to 


(2. 

0, 


A  ~  D  ±  R 
2C 

=  0 


A  =  ~l/K 


(26) 


of  which  the  first  is  truly  familiar  from  standard  cavity 
mode  analysis  while  the  last  is  from  dispersive  system 
theory. 

The  ray-pulse  matrix  theory  as  described  here  does  not 
include  active  modelockers  and  other  time-varying  ele¬ 
ments.  It  is  clear  that  the  theory  can  be  expanded  to  cover 
time  varying  elements  by  allowing  additional  elements  of 
the  matrix  to  vary.  Including  time-varying  modulation 


KOSTENBAUDER  RAY  PULSE  MATRICES  FOR  OPTICAL  SYSTEMS 


*  1157 


elements  should  make  it  possiMe  to  determine  the  true 
combined  temporal  and  spatial  pulse  shape  of  an  actively 
mode-locked  laser;  but  we  are  already  in  a  position  to  see 
that  the  correct  configuration  for  a  passively  modelocked 
system  will  have  A  near  0  as  this  means  tliat  the  pulse 
does  not  spread  much  during  its  propagation  around  the 
cavity,  and  thus  makes  easy  work  for  the  nonlinear  mech¬ 
anisms  that  lead  to  pulse  shortening. 

VIII.  Alternate  Ray-Pulse  Formulation 

For  compatibility  with  earlier  forms  of  ray  matrices,  we 
defined  the  ray-pulse  matrices  in  the  form  given  in  (4). 
An  alternative  approach  would  be  to  rearrange  the  vari¬ 
ables  with  both  coordinates  first,  both  frequencies  second 
in  the  form 


The  2x2  matrices  A,  B,  C,  D  are  then  obviously  biock- 
matrix  analogs  of  the  conventional  elements  A.  B.  C,  D. 

One  advantage  of  this  formulation  is  that  the  Gaussian 
beam  transformation  rule  can  be  written  in  the  simplified 
form 

Q,k..  =  IAQ,„  +  B|[CO,n  +  Dl'  (28) 

which  has  an  obvious  similarity  to  the  conventional  rule 
that  =  {Aq,„  +  B)/(Cq,„  +  D) . 

IX.  Conclusion 

We  have  shown  that  time-invariant  dispersive  optical 
systems  can  be  modeled  in  a  very  straightforward  way  by 
using  4x4  ray-pulse  matrices,  and  that  this  formulation 
allows  the  construction  of  a  Huygens  type  integral  oper¬ 
ator  that  describes  the  dispersive  system  up  to  quadratic- 


phases  in  both  spatial  and  temporal  frequencies.  The  abil¬ 
ity  of  4  X  4  matrices  to  convert  to  this  integral  operator 
form  relies  upon  the  fact  that  the  variables  come  in  con¬ 
jugate  pairs.  Several  extensions  of  this  work  are  possible. 
One  extension  would  be  to  include  elements  such  as 
Gaussian  apertures  whose  ray-pulse  matrices  would  have 
complex  entries.  Another  would  be  the  introduction  of 
time-varying  elements  such  as  modelockers.  By  including 
these  active  elements  it  should  be  possible  to  find  true 
eigenwaves  for  an  arbitrary  mode-locked  system,  with  re¬ 
sults  similar  to  the  Kuizenga-Siegman  |71  theory  of 
modelocking.  The  theory  could  also  be  extended  to  han¬ 
dle  off-axis  components,  and  finally,  cubic  and  higher  or¬ 
der  theories  might  be  developed  in  analogy  with  the  higher 
order  theories  of  monochromatic  systems.  For  the  disper¬ 
sive  problem,  there  is  an  interesting  quadratic  theory  as 
there  is  no  symmetry  which  rules  out  .second-order  terms 
as  there  is  in  the  familiar  case.  Also,  both  transverse  di¬ 
mensions  could  be  included  by  increasing  the  size  of  the 
matrix  to  6  x  6.  as  it  is  only  necessary  to  add  the  second 
slope  and  position  as  coordinates. 
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What  is  PARAXIA? 

Paraxia  is  a  famil-  Macintosh  applications  de¬ 
veloped  at  Stanfoni  urivcrsity  for  the  computa¬ 
tion  and  analysis  of  laser  resonator  modes  and  laser 
beam  propagation  problems.  Taking  advantage  of 
the  fast  processors  and  excellent  graphics  capabil¬ 
ities  of  modem  personal  computers,  together  with 
optimized  algorithms  iiKluding  Fast  Fourier  and  Fast 
Hankel  transforms  and  the  virtual  source  formalism, 
Paraxia  offers  a  set  of  powerful  analysis  and  design 
tools  that  provide  rapid  and  accurate  results  and  are 
extremely  easy  to  use. 

Paraxia  can  handle  in  general  any  separable 
paraxial  optical  system  in  cartesian  or  cylindrical  co¬ 
ordinates,  including  complex-valued  and  misaligned 
ray  matrices,  with  full  diffraction  effects  between 
apertures.  It  includes  the  programs; 

•  ABCD  provides  complex-valued  ray-matrix  and 
Hermite-gaussian  mode  analyses  for  arbitrary 
paraxial  resonators  and  optical  systems,  includ¬ 
ing  astigmatism  and  mis^ignment  in  each  ele¬ 
ment 

•  FRESNEL  uses  FFT  and  FHT  methods  to  propa¬ 
gate  an  arbitrary  wavefront  through  any  arbitrary 
paraxial  optical  system  using  Huygens’  integral 
in  rectangular  or  radial  coondinates.  The  wave- 
front  can  be  multiplied  by  an  arbitrary  mask 
or  mirror  profile  before  or  after  each  succes¬ 
sive  propagation  through  the  system;  a  scripting 
capability  makes  iterative  Fox  and  Li  resonator 
calculations  easy  to  implement. 

•  VSOURCE  implements  the  virtual  source  anal¬ 
ysis  to  calculate  accurate  eigenvalues,  eigen- 
modes,  excess  noise  factors  and  cross-power  co¬ 
efficients  for  unstable  resonators  with  both  cir¬ 
cular  mirrors  and  rectangular  hard-edged  mirrors 
(including  misaligned  rectangular  systems)  over 
a  full  range  of  equivalent  Fresnel  numbers. 

•  VRM  carries  out  mode  calculations  and  design 
analyses  for  gaussian  variable-reflcctivity-mirror 
lasers. 


How  to  get  Paraxia? 

Paraxia  is  available  from  Stanford  University  at: 
Software  Distribution  Center 
857  Serra  Street 
Stanford,  CA  94305-6225 
(415)  723-0651 


Why  use  PARAXIA? 

Because  PARAXIA  is  general,  fast,  accurate,  and  easy 
to  use. 


•  Wide  range  of  applications 

Paraxia  can  handle  arbitrary  paraxial  optical 
systems  with  only  a  few  restrictions,  and  can  be 
used  to  solve  a  wide  variety  of  problems,  from 
resonator  mode  calculations  to  far-field  beam 
patterns  to  misalignment  sensitivity  analysis. 

•  Performance 

Paraxia  uses  optimized  algorithms  for  fast  and 
accurate  results.  (FRESNEL  propagates  a  wave- 
front  of  1024  points  through  the  system  in  3 
seconds/iteration.  VSOURCE  computes  10  sets 
of  lowest  and  higher-order  unstable-resonator 
eigcnvalues/second.) 

•  User-friendiy  interface 

All  Paraxia  programs  are  highly  interactive 
using  an  efficient  user-friendly  interface  con¬ 
sistent  with  Macintosh  User  Interface  guide¬ 
lines.  Menus,  dialog  boxes  and  help  files  make 
Paraxia  easy  to  use  even  for  researchers  not 
familiar  with  the  Macintosh  computer. 

•  Compatibility 

Each  Paraxia  program  can  produce  its  own 
graphic  output  and  exchange  data  files  with  other 
Paraxia  programs;  but  all  the  programs  can  also 
export  data  files  to  a  variety  of  other  Macintosh 
applications,  including  major  graphics  packages, 
for  further  analysis  or  enhanced  display. 
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ABCD 


Ray-matrix  algebra  made  easy 


ABCD  is  a  highly  interactive  program  that 
allows  you  to  “build”  an  optical  system  as  a 
succession  of  optical  elements  (lenses,  mirrors, 
Brewster  plates,  etc.),  and  edit  it  in  the  same 
way  you  edit  text  on  the  Macintosh,  using  the 
Cut,  Copy,  Paste,  Clear  and  Undo  commands. 
The  various  parameters  (focal  length,  thickness, 
etc.)  specified  independently  for  each  element 
can  be  edited  at  any  time,  can  be  made  variable, 
and  each  element  can  be  made  misaligned  and/or 
astigmatic.  ABCD  will  propagate  an  arbitrary 
gaussian  beam  (e.g.,  the  eigenmode)  through  the 
system,  and  calculate  the  variation  of  the  gaus¬ 
sian  beam  parameters  (radius  of  curvature  and 
width)  at  a  given  location  when  varying  the  pa¬ 
rameters  of  the  optical  elements. 

Examples  of  problems  ABCD  can  easily 
solve  include: 

•  the  design/analysis  of  stable  resonators  (ei¬ 
genmode  calculation,  astigmatism  correc¬ 
tions,  regions  of  stability,  etc.) 

•  the  design/analysis  of  coupling  optics  (e.g., 
coupling  a  laser  beam  into  an  optical  fiber.) 

•  the  calculation  of  an  overall  ABCD  matrix 
to  be  used  in  FRESNEL. 

In  addition  to  being  a  powerful  investigation 
tool,  ABCD  also  proves  very  useful  as  a  didactic 
software  for  teaching  paraxial  optics. 

■  -■ 

- 1 

Each  element  Elements  can  be  astigmatic 

can  be  edited  at  any  time.  (like  this  ID  GRIN  slab). 


location)  when  varying  the  parameters 
of  the  optical  elements. 


A  second  window  shows  the  propagation 

of  an  arbitrary  gaussian  beam  through  the  system 

(both  transverse  directions  are  shown). 


Stanford  Resonator  arxt  Optics  Programs  —  FRESNEL 

General-purpose  wavefront  propagation 


FRESNEL  will  propagate  an  arbitrary  wave- 
front  successively  through  a  hard-edg^  aper¬ 
ture,  an  arbitrary  mask,  and  an  arbitray  paraxial 
system  described  by  a  complex-valued  ABCD- 
EF  matrix.  An  iterative  option  allows  for  Fox 
&  Li  calculations  of  resonator  eigenmodes. 

FRESNEL  uses  Fast  Fourier  and  Fast  Hankel 
transform  algorithms  for  fast,  accurate  results, 
even  with  large  numbers  of  sampling  points. 
FRESNEL  has  import  and  export  capabilities,  to 
exchange  wavefront  data  file  with  other  pro¬ 
grams,  and  can  be  run  either  in  interactive  or 
batch  mode  using  a  script  file  to  describe  longer 
jobs. 


Examples  of  problems  FRESNEL  can  easily 
solve  include: 

•  the  propagation  of  a  non-gaussian  wavefront 
through  an  arbitrary  paraxial  system. 

•  finding  the  eigenmodes  of  any  type  of  res¬ 
onator  (stable  or  unstable). 

The  hard-edged  aperture  and  arbitrary  mask 
in  front  of  the  system  allow  for  modeling  of  such 
elements  as: 

•  hard-edged  mirrors 

•  hole  coupling 

•  supergaussian  variable-reflectivity-mirrors 

•  gain  sheets 


The  input  wavefront  can  either  be 
a  build-in  analytic  function  or 
an  imported  user-defined  data  file. 


ie  Edit  Define  Run  illindouis 


Explicit  dialog  boxes 
allow  for  a  highly 
interactive  interface. 

1 


The  output  coordinate 
can  be  scaled  within 
the  calculation  process. 


On-line  help 
makes  the 
programs 
even  easier  to  use 


The  output  wavefront  can  be  exported 
directly  to  graphics  programs, 
such  as  Igor™. 


The  fast  FFT  and  FHT  algorithms 
allow  for  calculations  with  large 
numbers  of  sampling  points 
in  a  fairly  short  time. 

1024  points:  3  seconds/iteration 
16384  points:  50  seconds/iteration 


Stanford  Resonator  and  Optics  Programs 


^ograms  .VSOURCE 

All  about  unstable  resonators 


VSOURCE  implements  the  virtual  source 
formalism  for  the  study  of  hard-edged  strip  or 
circular  unstable  resonators.  VSOURCE  will 
compute  in  real  time  multiple  eigenvalues  as  a 
function  of  the  effective  Fresnel  number  and/or 
the  magnification  of  the  resonator,  as  well  as  the 
corresponding  eigenwaves  and  associated  over¬ 
lap  integrals  (excess  noise  factor,  cross-power, 
etc.). 


The  beam  profile  at  other  locations  inside 
the  resonator,  as  well  as  the  far-field  pattern 
outside  the  resonator,  can  be  calculated  by  ex¬ 
porting  the  data  to  FRESNEL.  In  cartesian  coor¬ 
dinates,  VSOURCE  can  handle  asymmetric  (or 
misaligned)  resonators,  allowing  the  study  of 
one-sided  unstable  resonators. 
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Eigenvalues  are  displayed 
both  in  modulus  and 
in  the  complex  plane. 


EKcess  Noise 


Eigenwaves  can  be 
exported  to  FRESNEL 
or  to  graphics  programs. 


Stanford  Resonator  and  Optics  Programs  _ VRM 

Gaussian  varlable-reflectivity-mirror  resonators 


VRM  is  a  numerical  implementation  of  complicated  analyt¬ 
ical  results  obtained  for  the  gaussian  variable-reflectivity-mirror 
resonators  (both  stable  and  unstable).  VRM  will  compute  the  vari¬ 
ation  of  the  eigenvalues  as  a  function  of  the  equivalent  (gaussian 
aperture)  Fresnel  number  and/or  the  magnification  of  the  resonator. 
For  selected  eigenvalues,  VRM  will  compute  the  corresponding 
eigenwaves  at  the  output  mirror  plane,  and  associated  overlap  in¬ 
tegrals  (excess  noise  factor,  cross-power,  etc.). 
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ABSTRACT 


There  is  a  continuing  interest  in  the  design  and  calculation  of  laser  resonators  and  optical  beam 
propagation.  In  particular,  recently,  interest  has  increased  in  developing  concepts  such  as  one-sided 
unstable  resonators,  supergaussian  reflectivity  profiles,  diode  laser  modes,  beam  quality  concepts, 
mode  competition,  excess  noise  factors,  and  nonlinear  Kerr  lenses. 

To  meet  these  calculation  needs,  I  developed  a  general-purpose  software  package  named 
PaRAXIA^,  aimed  at  providing  optical  scientists  and  engineers  with  a  set  of  powerful  design  and 
analysis  tools  that  provide  rapid  and  accurate  results  and  are  extremely  easy  to  use.  ParaxiA  can 
handle  separable  paraxial  optical  systems  in  cartesian  or  cylindrical  coordinates,  including  complex¬ 
valued  and  misaligned  ray  matrices,  with  full  diffraction  effects  between  apertures.  It  includes  the 
following  programs: 

ABCO  provides  complex-valued  ray-matrix  and  gaussian-mode  analyses  for  arbitrary  paraxial 
resonators  and  optical  systems,  including  astigmatism  and  misalignment  in  each  element.  This 
program  required  that  I  generalize  the  theory  of  gaussian  beam  propagation  through  real-valued  ray 
matrices  to  the  case  of  an  off-axis  gaussian  beam  propagating  through  a  misaligned,  complex-valued 
ray  matrix. 

FRESNEL  uses  FFT  and  FHT  methods  to  propagate  an  arbitrary  wavefront  through  an  arbitrary 
paraxial  optical  system  using  Huygens’  integral  in  rectangular  or  radial  coordinates.  The  wavefront 
can  be  multiplied  by  an  arbitrary  mirror  profile  and/or  saturable  gain  sheet  on  each  successive 
propagation  through  the  system.  Among  other  calculations,  I  used  FRESNEL  to  design  a  one-sided 
negative-branch  unstable  resonator  for  a  free-electron  laser,  and  to  show  how  a  variable  internal 
aperture  influences  the  mode  competition  and  beam  quality  in  a  stable  cavity. 

VSOURCE  implements  the  virtual  source  analysis  to  calculate  accurate  eigenvalues,  eigenmodes, 
excess  noise  factors  and  cross-power  coefficients  for  unstable  resonators  with  both  circular  and 
rectangular  hard-edged  mirrors  (including  misaligned  rectangular  systems).  I  used  VSOURCE  to 
show  the  validity  of  the  virtual  source  approach  (by  comparing  its  results  to  those  of  FRESNEL),  to 
study  the  properties  of  hard-edged  unstable  resonators,  amd  to  obtain  numerical  values  of  the  excess 
noise  factors  in  such  resonators. 

VRM  carries  out  mode  calculations  and  design  analyses  for  gaussian  variable-reflectivity-mirror 
lasers.  It  implements  complicated  analytical  results  that  I  derived  to  point  out  the  large  numerical 
value  of  the  excess  noise  factor  in  geometrically  unstable  resonators. 
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The  self-filtering  unstable  resonator  (SPUR)  can  provide  excellent 
mode  properties  for  high-gain  lasers  having  active-medium  Fresnel  num¬ 
bers  up  to  a  few  times  unity  or  perhaps  a  few  tens.  This  note  is  merely 
to  point  out  that  the  SPUR  design  is  generally  not  suitable  for  much 
larger  Fresnel  number  devices,  because  the  output  coupling  that  results 
for  single-mode  operation  with  full  energy  extraction  then  becomes  unac¬ 
ceptably  large. 
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The  self-filtering  unstable  resonator  (SPUR)  design  shown  in  Figure  1  has  been 
proposed  and  experimentally  demonstrated  by  Gobbi  and  Reali  [1-5]  and  a  number 
of  other  workers  [6-1  Ij  in  recent  years.  In  essence  the  SPUR  is  a  negative-branch 
confocal  unstable  resonator  in  which  a  small  aperture  through  the  scraper  mirror 
also  provides  an  internal  spatial  filter  at  the  common  focal  point  of  the  two  mirrors. 
As  this  earlier  work  hats  well  demonstrated,  the  SPUR  cavity  is  a  simple  and  practi¬ 
cal  design  which  can  provide  excellent  mode  properties,  including  very  good  beam 
quality  and  good  energy  extraction,  for  a  wide  variety  of  moderate  to  high-gain 
lasers.  It  is  therefore  a  very  practical  and  useful  design  for  many  laiser  systems. 

The  one  additional  point  that  I  would  like  to  make  in  this  comment  is  that 
despite  its  many  advantages  there  are  some  limitations  to  the  use  of  the  SPUR. 
This  type  of  resonator  is  associated  with  the  class  of  “unstable”  resonators  and 
might  therefore  seem  potentially  suitable  for  large  Fresnel  numbers  as  pointed  out 

\Z 

by  Anan’ev  [S].  What  I  wish  to  establish  in  this  comment  is  that  the  SPUR  is 
inherently  limited  to  laser  systems  with  at  best  moderate  values  of  the  leiser  Fresnel 
number,  that  is  to  Fresnel  numbers  on  the  order  of  a  few  times  imity,  up  to  a 
few  multiples  of  ten  at  the  most.  The  SPUR  design,  at  leeist  in  its  basic  form,  is 
not  suitable  for  laser  devices  in  which  the  inherent  Fresnel  numbers  are  much  larger 
than  a  few  multiples  of  ten  because  the  output  coupling  that  results  for  single-mode 
operation  with  full  energy  extraction  then  becomes  unacceptably  large. 

This  point  can  be  demonstrated  as  follows.  Figure  1  shows  the  basic  parameters 
of  the  SPUR  design  (the  notation  and  formulas  are  taken  from  the  paper  of  Boffa 
et  ai  [9]).  This  and  earlier  references  have  demonstrated  that  if  the  self-filtering 
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aperture  in  the  output  coupling  mirror  has  an  optimum  diameter  given  by 


d  =  2a«2  X  (0.61/2A)‘/2 


(1) 


then  this  resonator  will  produce  a  very  clean  near-gaussiein  lowest-order  mode  with 
a  gaussian  spot  size  w  «  Ma,  where  M  =  /i//2  is  the  geometric  magnification 
of  the  unstable  resonator.  The  magnitude  of  the  eigenvalue  7  for  the  lowest-order 
mode  and  thus  the  power  loss  per  round  trip  1  —  j-yp  will  then  be  given  to  a  good 
approximation  by 


1  -  jyp  =  1  -  il  w  1  - 


/ 1.968  \ 

0.6084' 

V  ) 

^  M2 

(2) 


The  reflectivity  value  R  =  (7P  can  be  viewed  as  the  effective  power  reflectivity  of 
the  output  mirror  or  output  coupler  in  this  cavity. 


An  additional  parameter  that  is  of  large  practical  importance  in  high-power  or 
high-energy  laser  designs  is  the  Fresnel  number  [H]  which  characterizes  the  laser 
medium  from  which  the  laser  mode  is  to  extract  power  or  energy.  In  particular,  if 
2A  is  the  diameter  and  L  the  length  of  the  laser  gain  medium,  then  it  is  convenient 
to  define  a  Fresnel  number  Nt  given  by 

(3) 

This  is  sometimes  called  the  “tube  Fresnel  number”  because  it  characterizes  the 
laser  tube  or  laser  rod  from  which  the  laser  energy  is  to  be  extracted.  If  one  is  to 
obtain  significant  energy  extraction  from  a  laser  gain  medium  of  diameter  2A  using 
a  gaussian  beam  profile  with  mode  spot  size  tv,  while  avoiding  significant  clipping 
of  the  gaussian  mode  by  the  tube  edges,  the  mode  size  and  the  tube  diameter  need 
to  be  related  by  a  ratio  on  the  order  of 


2A 


JTW 


(4) 
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For  the  paxticvilax  case  of  a  SFUR  cavity,  the  tube  Fresnel  number  Nt  and  the 
magnification  M  will  then  be  related  by 

2  2 

iVr«(f)  1.5  X  (/,/£)  X  M  (6) 

Since  Nt  is  determined  by  the  laser  medium,  and  M  is  more  or  less  constrained 
by  the  allowable  output  coupling  through  Eq.  (2),  the  only  adjustable  parameter  is 
the  ratio  of  the  mirror  focal  length  /i  (which  is  also  to  a  first  approximation  the 
overall  resonator  length)  to  the  length  L  of  the  gziin  medium. 

The  most  difficult  situations  for  laser  resonator  design  arise  in  high-power  lasers 
having  Izirge-diameter  gain  media  and  therefore  values  of  Nt  that  are  sometimes 
very  much  greater  than  unity.  A  laser  tube  with  bore  diaimeter  2A  =  1  cm  and 
length  £  =  1  m  at  a  visible  wavelength  A  =  500  nm,  for  example,  has  Nt  ^  50. 
Equation  (5)  indicates  that  efficient  single-mode  energy  extraction  from  such  large- 
Nt  laser  media  can  only  be  accomplished  with  an  SFUR  design  either  by  using 
a  very  large  magnification  M,  which  means  very  large  output  coupling  and  hence 
requires  very  high  Itiser  gain,  or  by  using  a  mirror  focal  length  fi  very  much  longer 
than  the  gain  medium  length  £,  or  possibly  by  using  a  folded  multipath  axrangment 
through  the  laser  medium.  The  latter  arrangments  tend  to  be  both  mechanically 
^lnd  optically  inconvenient. 

To  demonstrate  this  requirement  in  a  graphic  fashion,  Figure  2  plots  the  re¬ 
lationship  between  the  tube  Fresnel  number  Nt  and  the  effective  output  mirror 
reflectivity  R  =  j-yp  for  three  different  ratios  of  the  mirror  focal  length  /i  to  the 
length  L  of  the  active  leiser  medium,  assuming  a  SFUR  design  which  satisfies  both 
the  single- mode  condition  of  Eq.  ( 1 )  and  the  full  power  extraction  condition  of  Eq. 
(4).  (Note  that  the  magnification  M  is  a  variable  parameter  along  each  of  these 
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curves.)  Note  that  in  each  case  the  effective  reflectivity  of  the  output  mirror  drops 
rapidly  to  values  below  10%  as  the  value  of  Nt  increases;  and  that  only  a  few  lasers 
operate  well  with  output  reflectivities  of  i?  =  10%  or  lower.  Most  of  the  high-power 
and  high-energy  lasers  of  practical  interest  operate  best  with  output  couplings  of 
20%  to  50%  or  larger  (though  there  axe  of  course  some  exceptions  to  this). 

The  general  observation  from  this  figure  is  that  the  SPUR  can  provide  rezison- 
able  output  coupling  values  or  mirror  reflectivities  for  leiser  gain  media  characterized 
by  values  of  the  Fresnel  number  Nt  ranging  from  a  few  times  unity  up  to  perhaps 
Nt  =  10  or  20.  The  SPUR  approach  cannot,  however,  provide  a  solution  for  single¬ 
mode  operation  with  good  energy  extraction  in  the  more  difficult  case  of  gain  media 
with  much  larger  tube  Fresnel  numbers,  because  the  required  values  of  magnifica¬ 
tion  needed  to  fill  the  gain  medium  will  make  the  effective  output  coupling  from 
the  cavity  lau'ger  than  can  be  tolerated.  Lasers  with  larger  values  of  Nt  will  require 
the  use  of  some  combination  of  either  the  original  hard-edged  unstable  resonator 
approach  [13]  or  some  form  of  variable-reflectivity-mirror  (VRM)  resonator  design 
[5,14]  or  possible  some  form  of  injection  locking  or  seeding  [7],  along  with  possible 
folding  or  multipass  operation  of  the  laser  cavity. 
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FIGURE  CAPTIONS 


Fig.  1.  Beisic  self-filtering  unstable  resonator  (SFUR)  design.  The  gain  medium  is 
assumed  to  have  an  outer  diameter  2 A  and  to  extend  over  a  length  L  <  fi. 

Fig.  2.  Effective  power  reflectivity  R  =  I7P  of  the  output  coupling  element  versus 
Fresnel  number  Nt  =  A^fLX  for  SFUR  designs  with  three  different  values  of 
mirror  focal  length  /i  to  active  gain  medium  length  L.  (The  magnification  M 
is  a  variable  parameter  along  each  of  these  curves. ) 
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In  the  retrospective  part  of  this  talk  I’d  like  to  take  a  look  backward  at  the 
many  important  and  interesting  concepts  in  optical  resonators  that  have  evolved 
since  lasers  were  first  developed  some  three  decades  ago. 

Then  in  the  prospective  part  I’ll  take  a  look  forward  at  some  of  newer  ideas 
that  are  now  emerging,  or  that  have  recently  been  revitalized,  and  that  are  likely 
to  play  an  important  role  in  the  coming  decades. 

Except  in  few  cases  I  will  not  try  to  identify  the  people  who  made  all  these 
contributions.  Let  me  apologize  in  advance  to  those  whom  I  do  not  mention,  and 
also  to  those  whose  favorite  ideas  I  may  have  overlooked. 

But  I  would  particularly  like  to  Jacques  Amaud  who,  along  with  several 
others,  made  helpful  suggestions  for  this  talk;  the  Air  Force  Ofiice  of  Scientific 
Research  which  has  supported  nearly  all  of  my  research  in  this  area  over  the 
years;  and  the  many  students  and  other  colleagues  with  whom  it  has  been  such 
a  pleasure  to  work  on  these  resonator  ideas  over  the  past  several  decades. 
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An  overall  list  of  the  most  important  ideas  that  have  been  developed  in  optical 
resonator  theory  and  practice  would  surely  include  all  of  the  following  areas: 


LASER  RESONATORS:  RETROSPECTIVE 

•  The  existence  of  modes  in  open-sided  resonatore 

•  Ray  matrices  and  gaussian  beams 

•  Gaussian  resonator  modes 

•  Generalized  paraxial  wave  optics 

•  Unstable  optical  resonators 

•  Biorthogonality  properties  of  optical  resonators 

•  Frequency-selective  and  tunable  laser  cavities 

•  Other  beam  and  resonator  ideas 


We’ll  try  to  take  a  more  detailed  look  at  some,  though  not  all,  of  these  topics. 
Let’s  walk  through  them  quickly,  one  by  one. 
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•  • 

Some  of  the  important  basic  ideas  that  emerged  in  the  first  days  of  the  laser 
era  included  the  following; 

Existence  of  “Modes”  in  Open-Sided  Resonators 

o  The  Fabry-Perot  laser  cavity  concept  (Schawlow  and  Townes) 
o  The  first  Fox  and  Li  recirculating-pulse  calculations 
0  Gobau’s  beam  waveguide  ideas 
o  The  equivalent  periodic  lensguide  viewpoint 
o  Existence  of  unique  resonator  eigenmodes  and  eigenvalues 
o  Existence  of  low-loss  curved-mirror  resonators 
o  Importance  of  the  resonator  Fresnel  number 
o  Transverse  modes  and  transverse  mode  beats 


The  ideas  of  ray  matrices  or  ABCD  matrices,  and  of  gaussian  beams  and 
resonators,  and  of  the  interconnections  between  them,  have  been  evolving  in 
tandem  ever  since: 


Ray  Matrices  and  Gaussian  Beams 

o  Ray  matrices  (ABCD  matrices) 
o  Periodic  focusing  systems  and  periodic  lensguides 
o  Ray  matrix  eigenvalues  and  eigenrays 
o  Stable  and  unstable  periodic  focusing  systems 
o  Hermite-gaussian  optical  beams 
o  Gaussian  beam  propagation  rules 

o  The  complex  radius  of  curvature  or  gaussian  q  parameter 

o  Imaginary  source  points,  the  confocal  parameter,  and  the  Rayleigh 
range 

o  Focal  phase  shifts  and  the  Gouy  effect 
0  Gaussian  ducts  and  graded-index  (GRIN)  rods  and  lenses 
o  Gaussian  mode  matching  and  the  Collins  chart 
o  Ray  matrices  and  Huygens*  integral 


Developing  the  analytical  and  experimental  imderstanding  of  gaussian  cav¬ 
ity  modes  and  their  application  in  lasers  has  been  an  important  part  of  this 
subject — and  also  one  of  the  many  areas  in  which  R.  V.  Pole  made  fundamental 
contributions: 


Gaussian  Resonator  Modes 

o  Stable  and  unstable  paraxial  optical  resonators 
o  Hermite-gaussian  and  Laguerre-gaussian  resonator  modes 
o  Relationship  to  periodic  focusing  systems 
0  The  resonator  stability  diagram 

0  The  confocal  cavity  concept 

-  Wieder  and  Pole’s  conjugate-concentric  laser  cavity 

-  Scanning  Fabry-Perot  interferometers 
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In  fact,  we’ve  now  come  to  vmderstand  that  not  just  geomerical  optics  but 
essentially  all  of  diffraction  optics  in  the  paraxial  approximation  can  be  incorpo¬ 
rated  into  a  generalized  complex-valued  ray  matrix  and  paraxial  optical  formu¬ 
lation  of  optical  beams  and  resonators: 


Generalized  Paraxial  Wave  Optics 

o  Gaussian  apertures  as  lenses  of  imaginary  focal  length 
o  Complex-valued  optical  elements  and  ABCD  matrices 
0  Complex-valued  paraxial  optical  systems 
o  Hermite-gaussian  functions  of  complex  argument 
0  Complex  spot  sizes  as  well  as  complex  q  parameters 

o  Generalized  paraxial  resonator  analysis 

-  Confined  and  unconfined  gaussian  modes 

-  Geometricaliy  stable  and  unstable  gaussian  modes 

-  Perturbation-stable  and  unstable  gaussian  modes 

-  Huygens  integral  and  complex-valued  ABCD  matrices 
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<  • 

The  tinstable  resonator  of  course  represented  a  different,  complicated,  chal¬ 
lenging,  and  very  useful  area  of  resonator  theory  and  practice: 


Unstable  Optical  Resonators 

o  Stable  modes  in  *^unstable”  systems. . .? 
o  Geometrical  mode  solutions  for  unstable  resonators 
o  Magnifying  and  demagnifying  spherical  wave  solutions 
o  Positive-branch  and  negative-branch  resonators 
o  Real  unstable  resonator  eigenmodes  and  eigenvalues 
o  FFTy  FHT  and  Prony  computational  methods 
o  Mode  crossings  and  the  equivalent  Fresnel  number 
o  Confocal  unstable  resonators 
o  Self-imaging  unstable  resonators 
o  Edge  waves,  and  the  Poisson  Spot  or  Spot  of  Arago 
o  The  asymptotic  approach,  or  ‘^virtual  source’*  analysis 
o  Self-focusing  (SFUR)  unstable  resonators 


Some  of  the  most  interesting  aspects  of  optical  resonators  arise  from  the  fact 
that,  unlike  most  physical  systems  with  whi^  we  are  familiar,  laser  cavities  axe 
generally  not  described  by  a  Hermitian  operator,  and  therefore  their  reson2Lnt 
modes  are  not  orthogonal — especially  the  modes  in  unstable  optical  resonators: 


Biorthogonality  Properties  of  Optical  Resonators 

o  Optical  resonator  modes  are  NOT  in  general  power  orthogonal 
0  More  power  in  one  mode  than  total  power  in  all  modes. . .? 
o  Matched  coupling  and  adjoint  coupling  methods 
o  Excess  quantum  noise  in  unstable  resonator  modes 
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There  are  many  other  intriguing  and  ingenious  ideas  related  to  laser  res¬ 
onators  that  have  sJso  been  developed  over  the  years,  in  order  to  solve  one  or 
another  of  the  practical  problems  facing  the  laser  designer: 


Frequency-Selective  and  Tunable  Laser  Cavities 
o  Multi- mirror  cavities 

o  Michelson,  Fox-Smith,  and  vernier  interferometer  cavities 
o  Tunable  birefringent  filters 
o  Prism  beam  expanders 
o  Etalon  mirrors 

o  Sagnac  interferometers  and  antiresonant  ring  cavities 
o  “Twisted-ribbon”  laser  cavities 
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And  still  more  clever  ideas,  inventions  and  concepts  that  have  been  applied 
to  less  conventionzd  sorts  of  laser  resonators  and  laser  beams: 


Other  Important  Beam  and  Resonator  Ideas 

o  Ring  laser  cavities  and  optical  diodes  and  isolators 
o  Optical  diodes  and  isolators 
o  Distributed-feedback  (DFB)  cavity  concepts 
o  Distributed  Bragg  reflector(DBR)  cavity  mirrors 
o  Waveguide  lasers  and  waveguide  laser  cavities 
o  Heterojunction  diode  laser  structures 
o  Corner-cube  and  rooftop  resonators 
o  Axiconsy  waxicons,  reflaxicons,  HSURIAs,  and  so-ons 
o  Whispering  gallery  modes  and  ^‘droplet**  resonators 
o  Pulse  compressors,  grating  and  prism  dispersive  systems 
o  Image  relay  systems  and  supergaussian  beam  profiles 
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Just  making  up  the  preceding  lists  of  topics,  without  examining  them  in 
detail,  shows  how  fm  we’ve  come,  and  how  mzmy  new  and  interesting  ideas  have 
emerged  in  the  resonator  field.  And  yet  the  topics  mentioned  so  fax  do  not  even 
include  all  of  the  complex  aspects  of  laser  oscillation  and  laser  mode  dynamics, 
such  as  are  illustrated  by  the  following  topics: 


Laser  Oscillation  and  Laser  Dynamics 

o  Laser  gain  saturation  effects 
o  Mode  build-up  and  mode  competition 
o  Spiking  and  relaxation  oscillations 
o  Laser  Q-switching 

o  Axial  mode  coupling  and  mode-locking 
o  Instabilities  and  chaotic  behavior  of  laser  oscillators 
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Nor  do  these  lists  include  the  even  more  complex  world  of  nonlinear  optical 
propagation  eflFectc,  and  optical  mode  behavior,  as  briefly  summarized  in  the 
following  list: 


Nonlinear  optical  propagation  effects 

o  Concept  of  phase  matching  in  nonlinear  optics 
o  Self-focusing  and  trapping 
o  Seif-phase-modulation 
o  Optical  bistability 
o  Optical  phase  conjugation 

o  Optical  chaos  in  passive  beams  and  resonators 
o  Solitons 
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But  now,  in  the  remaining  time,  let’s  turn  to  the  “prospective”  portion  of 
this  talk,  and  take  a  look  forward  at  some  newer  ideas  now  emerging  that  will 
play  an  important  role  in  coming  decades. 

Of  course,  some  of  these  “new”  ideas  are  in  fact  older  ideas  that  are  only 
now  becoming  really  important  in  practical  laser  designs,  while  others  represent 
ingenious  new  ideas,  developed  to  solve  new  problems. 

In  any  event,  I  beheve  a  broad  overview  of  the  ideas  now  emerging  in  the 
laser  resonator  field  would  surely  include  emerging  developments  in  each  of  the 
following  three  eu-eas: 


LASER  RESONATORS:  PROSPECTIVE 

•  New  kinds  of  resonators  for  new  kinds  of  lasers 

•  New  developments  in  diode  laser  resonators 

•  New  design  tools  and  beam  diagnostic  methods 


Again,  let’s  take  a  look  at  each  of  these  topics  in  more  detail. 


There  are  several  new  kinds  of  resonators  that  are  now  emerging  for  use  in 
new  kinds  of  laser  devices  that  have  great  promise  for  the  future.  A  list  of  these 
new  resonator  concepts  might  include  the  following,  and  I’ll  try  to  say  a  few 
words  about  each: 


New  Kinds  of  Resonators  for  New  Kinds  of  Lasers 

•  Tapered  reflectivity  and  VRM  laser  cavities 

-  Advances  in  tapered-reflectivity  mirror  fabrication 

-  Supergaussian  mirrors  and  mode  profiles 

•  Fibei>optic  resonators 

-  Couplers,  rings,  interferometers,  fiber  gyros 

•  Miniature  and  monolithic  laser  cavities 

-  Diode-pumped  solid-state  lasers 

-  Integrated-optical  resonator  structures 

-  “Platelet”  laser  structures 

•  Nonplanar  ring  resonators 

-  Monolithic  diode-pumped  NPR  lasers 

-  “UR-90”  single-sided  unstable  resonators 

•  Advances  in  high-power  slab  laser  designs 

•  Coupled  lasers  and  laser  oscillator  arrays 

•  Advances  in  FEL  resonator  designs 
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Of  course  one  of  the  most  important  and  active  areas  of  current  laser  research 
is  in  the  burgeoning  world  of  semiconductor  diode  lasers,  where  many  fascinating 
new  resonator  concepts  are  emerging  to  meet  the  special  needs  and  constraints 
of  this  class  of  leisers: 


New  Developments  in  Diode  Laser  Resonators 

•  Coupled-stripe  diode  laser  arrays 

•  Antiguidinge  (ROW)  diode  laser  stripe  arrays 

•  Talbot-plane  and  Talbot-mirror  resonator  concepts 

•  Grating-coupled  surface-emitting  diode  lasers 

•  Active  on-chip  waveguide  structures  and  MOPA  chains 

•  Stable-unstable  resonator  concepts  for  diode  lasers 

•  Vertical-cavity  surface-emitting  diode  laser  structures 

•  Etched-facet  diode  laser  mirrors 


Some  of  the  ideas  listed  above  will  turn  out  to  be  of  great  practical  importance 
in  the  future,  and  I  expect  some  of  them  will  be  only  of  historic^Ll  interest  a  few 
years  from  now.  We’ll  just  have  to  wait  and  see,  however,  which  of  these  ideeis 
fall  into  which  of  these  categories! 
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Finally,  I’d  like  to  say  a  few  words  about  three  current  developments  in  design 
tools  and  diagnostic  methods  related  to  laser  resonators  and  beams  in  which  I’ve 
been  fortunate  to  have  had  some  involvement: 


New  Design  Tools  and  Diagnostic  Methods 


•  New  measures  for  specifying  and  determining  laser  beam  quality 


•  Resonator  analysis  and  design  using  personal  computers 


•  Ray-pulse  matrices  for  dispersive  optical  systems 
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First  of  all,  after  talking  for  several  decades  about  “laser  beam  quality”  we 
are  finally  getting  Jiround  to  developing  real  instruments  that  will  really  measure 
real  laser  beam  quality.  In  fact  an  instrument  that  does  just  this  in  real  time  will 
be  on  display  for  the  first  time  at  this  meeting.  (It’s  been  developed  by  Coherent, 
Inc.,  if  I  may  be  allowed  mention  of  a  commercial  product.) 

At  the  same  time  we’re  developing  new  ways  of  defining  the  quality  of  an 
arbitrary  laser  beam,  and  making  theoretical  use  of  this  definition.  Part  of  this 
is  based  on  still  another  extension  of  ABCD  matrix  theory:  new  understanding  of 
how  the  moments  of  an  arbitrary  (nongaussian)  beam  can  be  propagated  through 
an  arbitrary  paraxial  system,  using  only  the  ABCD  matrix  elements. 


I’d  also  like  to  mention  briefiy  the  development  by  my  students  and  my¬ 
self  of  resonator  and  optical  beam  analysis  and  design  programs  that  can  be 
implemented  on  desktop  personal  computers  rather  than  mainframes.  We’ve 
developed  a  set  of  four  useful  programs  for  doing  laser  beam  propagation  and 
Fox-and-Li  type  resonator  calculations,  unstable  resonator  design  and  analysis, 
and  the  analysis  of  complex  ABCD  matrix  systems. 

These  programs  will  become  available  from  Stanford  University  some  time 
later  this  summer.  They  rim  on  Macintosh  computers,  and  can  do  nearly  every¬ 
thing  that  used  to  be  done  a  decade  or  more  ago  on  mainframe  supercomputers. 


Finally  I’d  like  to  mention  a  very  general  and  powerful  “ray-pulse  matrix” 
formulation  recently  developed  by  one  of  my  students.  Ad  Kostenbauder.  This 
formulation  extends  the  familiar  ray-matrix  formalism  to  include  both  spatial 
and  temporal  (or  spectral)  variations,  making  it  possible  to  handle  both  the 
time  and  the  space  dependence  of  optical  beams  and  pulses  in  one  calculation. 

This  requires  4x4  rather  th?m  just  2x2  matrices,  of  course,  but  this  approach 
has  already  proven  extremely  valuable  for  calculations  like  the  propagation  of 
femtosecond  pulses  in  dispersive  optical  systems  and  the  behavior  of  dispersive 
prism  and  grating  beam  expanders. 
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In  conclusion,  I  think  one  can  only  be  impressed  first  by  the  large  number 
of  clever  and  innovative  ideas  that  have  been  developed  in  this  field  over  the 
years,  and  second  by  the  often  remarkable  simplicity  and  elegance  of  the  basic 
underlying  principles  that  have  emerged.  It  really  seems  that  the  most  important 
ideas  are  often  the  simple  ones,  and  that  nature  rewards  us  with  elegance  in  the 
most  basic  ideas. 

I  think  R.  V.  Pole  recognized  these  principles,  and  practiced  them  in  his  own 
work,  and  I’m  honored  to  speak  today  in  his  memory. 

Looking  at  the  present  and  the  future,  one  can  also  still  find  great  vitality 
in  this  field,  with  a  set  of  challenging  new  problems  and,  one  hopes,  still  more 
interesting  and  elegeint  ideas  yet  to  come. 

I’d  like  to  thank  once  again  the  CLEO/IQEC  organizers  for  inviting  me  to 
speak,  especially  Steve  Guch  who  suggested  the  title  for  this  talk;  0Uid  the  many 
students  and  colleagues  who’ve  stimulated  me  and  contributed  to  my  work  over 
the  yezu^. 

And,  I’d  also  like  to  acknowledge  once  again  the  research  support  from  the 
General  Physics  division  of  the  Air  Force  Ofiice  of  Scientific  Research  (AFOSR), 
which  has  made  possible  almost  all  of  my  own  work  and  that  of  many  others  in 
this  field. 


18 


1146 


f 

IEEE  JOURNAL  OF  QUANTUM  ELECTRONICS.  VOL.  27.  NO  5.  MAY  1991 

Defining  the  Effective  Radius  of  Curvature  for  a  Nonideal  Optical  Beam 

A.  E.  Siegman 


Abstract— A  simple  expression  is  given  for  evaluating  the  ef¬ 
fective  or  average  radius  of  curvature  /Hz)  of  an  arbitrarily 
wrinkled,  distorted,  or  nondiffraction-limited  optical  beam. 
This  real-beam  radius  of  curvature  obeys  the  same  free-space 
propagation  equation  as  does  the  wavefront  curvature  of  an 
ideal  Gaussian  beam;  and  extracting  this  “correctable  curva¬ 
ture”  from  the  beam  at  any  plane  simultaneously  collimates  the 
beam,  makes  that  location  a  waist  for  the  real  beam,  and  min¬ 
imizes  the  far-held  beam  divergence  starting  from  that  posi¬ 
tion. 


Fig.  I.  An  arbitrary  nondilTraction-limited  optical  beam  propagating  out¬ 
ward  from  an  input  plane  c  =  c,  towards  a  waist  at  c  =  c„. 


CONSIDER  a  monochromatic  optic  beam  propagating 
within  the  paraxial  approximation,  but  suppose  that 
this  beam  is,  in  general,  neither  spherical  nor  collimated 
nor  diffraction-limited.  To  say  this  in  another  way,  the 
wavefront  of  the  beam  across  any  transverse  plane  may 
be  wrinkled  as  well  as  curved,  and  the  transverse  ampli¬ 
tude  profile  may  have  complicated  and  irregular  varia¬ 
tions  as  well.  What,  then,  is  an  appropriate  definition  for 
the  effective  or  averaged  or  “removable”  radius  of  cur¬ 
vature  R  (z)  of  this  real-beam  wavefront  at  that  transverse 
plane?  The  purpose  of  this  note  is  to  show  hpw  to  define 
an  effective  radius  of  curvature  R(z)  for  such  a  real  or 
nonideal  laser  beam. 

By  way  of  analytical  background,  let  the  complex  pha- 
sor  amplitude  profile  £(x,  z)  of  this  beam  be  specified  at 
a  plane  z  =  Zt  which  is  not  necessarily  coincident  with 
the  waist  location  z  =  Zo  of  the  beam,  as  illustrated  in 
Fig.  1.  The  phasor  amplitude  £{x,  z)  of  the  beam  at  any 
plane  z  can  then  be  related  to  a  complex  phasor  amplitude 
P(s,  z)  which  describes  the  beam  in  the  angular  or  spatial- 
frequency  domain  through  the  Fourier  transform  relation 

Joo 

Pis,  z)e~^^"' ds  (I) 

-00 

with  Pis,  z)  given  by  the  reverse  transform 

JOO 

£ix,  dx.  (2) 

—  00 

The  spatial  frequency  variable  s  is  related  to  the  angular 
direction  0  of  a  corresponding  plane-wave  component  of 
the  beam  by  5  ®  '  sin  fl,  which  becomes  $  ~  Xs  in  the 

paraxial  approximation.  The  spatial-frequency  distribu¬ 
tion  Pis.  z)  then  propagates  in  free  space  according  to 

Pis,  z)  =  Pis,  (3) 


which  shows  that  individual  components  of  Pis,  z)  rotate 
in  phase,  but  do  not  change  in  magnitude  with  distance  z- 
We  also  assume  for  simplicity  that  £ix,  z)  is  normalized 
so  that  j“a,  l£ix,  z)(^  dx  =  j”*  ds  =  1. 

It  can  then  be  shown  that  the  variance  a^iz)  of  this  real- 
beam  profile  in  the  transverse  direction,  evaluated  using 
the  intensity  distribution  |£(x,  z)|^,  will  have  an  axial 
variation  given  by 

<j^(z)  =  a^(Zi)  -  /4i  X  (z  -  Zi)  +  X  (z  -  Z|)^ 

(4) 

where  a^izi)  is  the  spatial  variance  at  the  input  plane  z  = 
Zi;  ffj  is  the  variance  of  the  beam  in  the  angular  or  spatial- 
frequency  domain,  evaluated  using  the  spatial-frequency 
intensity  distribution  |^(s,  z)p;  and  /4|  is  the  value  of  the 
function 


evaluated  at  z  =  Zj.  The  Fourier  transform  relations  be¬ 
tween  £ix,  z)  and  Pis,  z)  imply  that  there  will  be  corre¬ 
sponding  Fourier  transform  relations  between  the  quan¬ 
tities  d£ix,  z)/dx  and  -jlxsPis,  z),  and  also  between 
the  quantities  xEU,  z)  and  -ij/2x)  dPis,z)/ds.  If  these 
Fourier  transform  relations  are  substituted  into  (5),  the 
function  A  (z)  can  also  be  written  as 


X  f°°  r 


d£*ix.  z) 
dx 


-  £*ix.  z) 


d£ix,  z) 
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with  A,  again  being  this  function  evaluated  at  z  =  z,. 

The  axial  variation  of  a‘(z)  in  (4)  can  be  rewritten  as 
an  expansion  about  a  waist  location  z  =  Zo  in  the  form 

flj.  [2J 
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a  (^)  =  a5  X  1  + 


where  al  is  the  minimum  value  of  the  spatial  variance  at 
the  real-beam  waist,  given  by  al  =  a'(zi)  - 
and  Zo  is  the  waist  location  which  is  related  to  by  zo  = 
Zi  -I-  /4|/2X'aj.  It  is  also  convenient  to  introduce  a  beam 
quality  factor  {2),  [3]  which  is  related  to  the  space-beam- 
width  product  aoUj  for  the  real  beam  by  M‘  =  4iraoa,. 
This  quality  factor  has  the  property  that  >  1  for  any 
real  beam,  reducing  to  the  limiting  value  =  1  only 
for  an  ideal  TEMqo  Gaussian  beam.  Equating  the  quad¬ 
ratic  terms  in  (3)  and  (7)  then  gives  a  formula  for  the  Ray¬ 
leigh  range  Zr  of  the  real  beam,  namely. 


This  is  the  same  as  the  formula  for  the  Rayleigh  range  of 
an  ideal  TEMqq  Gaussian  beam  with  waist  variance  oq, 
which  is  Zr  =  TTwl/X  =  4iral/\,  except  that  the  Rayleigh 
range  Zr  for  a  real  nonideal  beam  is  reduced  by  the  factor 
>  1.  Equation  (7)  confirms  that  any  arbitrary  real 
beam  always  has  a  waist,  with  a  spot  size  or  variance 
which  varies  quadratically  on  either  side  of  the  waist  just 
like  an  ideal  Gaussian  beam,  except  that  the  Rayleigh 
range  for  the  real  beam  is  reduced  by  the  beam  quality 
factor 

A  general  expression  for  the  effective  radius  of  curva¬ 
ture  R  (z)  of  any  real-beam  wavefront  at  any  arbitrary  plane 
z  can  then  be  written  as 

—  -  ^  (9) 

R(Z)  2a\z) 

where  A(z)  can  be  evaluated  using  either  (5)  or  (6).  This 
expression  may  not  be  entirely  new.  Essentially  this  same 
formula  was  used  on  an  empirical  basis  as  a  way  of  ex¬ 
tracting  the  macroscopic  wavefront  curvature  in  early  nu¬ 
merical  beam  propagation  calculations  [4],  and  the  same 
formula  has  also  arisen  in  recent  work  by  Belanger  [5]. 
The  primary  objective  of  this  note,  in  light  of  current  dis¬ 
cussions  of  laser  beam  quality  and  laser  beam  laser  beam 
propagation,  is  to  reiterate  the  usefulness  of  (9)  as  the 
basic  definition  of  effective  wavefront  curvature,  and  to 
point  out  three  or  four  different  ways  in  which  this  defi¬ 
nition  is  physically  meaningful  as  an  effective  radius  of 
curvature. 

The  physical  significance  of  the  definition  of  R(z)  given 
in  (9)  can  be  demonstrated  first  of  all  by  noting  that  this 
definition,  together  with  the  other  equations  above,  leads 
to  the  free-space  propagation  law 


R(z)  =  (z  —  Zo)  + 


(z  ~  Zo) 


where  Zr  is  the  same  real-beam  Rayleigh  range  as  in  (9). 
In  other  words,  with  the  use  of  the  real-beam  Rayleigh 
range,  both  the  real-beam  variance  a“(z)  from  (7)  and  the 
real-beam  effective  radius  R(z)  from  (9)  obey  exactly  the 


same  free-space  propagation  rules  as  do  the  correspond¬ 
ing  quantities  for  an  ideal  Gaussian  beam. 

Another  way  of  arriving  at  the  curvature  definition  in 
(9)  is  to  ask  what  is  the  radius  of  curvature  R(z)  for  a  real- 
beam  wavefront  at  any  plane  z  which,  when  extracted  from 
the  input  phasor  amplitude  £(jr,  z)  at  this  plane,  will  make 
the  modified  wavefront  at  this  plane  be  a  waist,  which 
means,  in  turn,  that  the  modified  value  of  A(z)  after  this 
extraction  should  be  zero.  This  would  seem  to  be  a  mean¬ 
ingful  way  of  defining  an  effective  radius  of  curvature  for 
the  real  beam  since  it  determines  what  radius  of  curvature 
has  to  be  extracted  from  the  real  wave  using  an  appropri¬ 
ate  thin  lens  to  convert  it  into  a  collimated  beam  waist. 
To  determine  this  removable  curvature,  we  can  suppose 
that  the  input  wave  E(x,  z)  at  plane  z  is  passed  through  a 
thin  lens  of  focal  length  /  so  that  the  beam  wavefront  is 
transformed  to 


E'(x,  z)  =  E(x,  z)  X  exp  +7—  (11) 

J  ^ 

where  £'(x^  z)  is  the  collimated  (although  potentially  still 
wrinkled)  phasor  amplitude  for  the  wave  that  will  be  left 
after  the  effective  curvature  has  been  removed  using  the 
lens  of  focal  length  /.  This  is  sometimes  referred  to  as 
removing  the  “correctable  divergence”  of  the  beam  at  the 
plane  z.  The  modified  value  of  the  integral  A  (z)  after  this 
curvature-extraction  process  will  be  given  by 

/i'(2)  =  -j^[°  x\£'^-^-e'*~]dx 


=  A(z)  + 


2ct'(z) 


Requiring  that  /!  '(z)  =  0  after  the  lens  then  gives 


/  /?(z)  2a'(z) 

exactly  as  expected. 

As  another  approach  to  the  definition  of  effective  cur¬ 
vature,  we  can  note  from  (4)  that  the  far-field  angular 
spread  of  an  arbitrary  beam  is  directly  proportional  to  the 
spatial-frequency  standard  deviation  a,.  The  effective  ra¬ 
dius  of  curvature  for  an  arbitrary  wavefront  would  then 
seem  to  be  that  radius  which,  if  extracted  from  the  wave- 
front  using  a  thin  lens  as  above,  will  minimize  the  result¬ 
ing  value  of  a;.  From  standard  Fourier  transform  identi¬ 
ties,  the  value  of  ctJ  can  be  calculated  using  either  of  the 
expressions 


5‘|  P{s,  Z)p  ds 


d£{x.  z) 
dx 


Let  the  complex  wavefront  E(.r,  z)  at  plane  z  again  be 
pas.sed  through  a  lens  of  focal  length  /  so  that  the  wave- 
front  is  modified  as  in  (1 1).  Putting  this  modified  wave- 
front  into  the  last  integral  in  (14)  shows  that  a;  after  the 
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The  value  of /which  minimizes  a',  is  then  exactly  the  same 
as  in  (13). 

Finally,  we  can  note  that  if  the  complex  phasor  ampli¬ 
tude  of  the  beam  is  separated  into  amplitude  and  phase 
functions  in  the  form 


£(x,  z)  *  |£(x,  z)|  exp  z)]  (16) 


then  the  general  definition  of  (9)  can  be  written  as 


1 

R{z) 


|£(jr,  z)p  dx. 


(17) 


This  says  that  I /R(z)  represents  a  weighted  average  of 
the  wavefront  curvature  term  [jc  3<(>(x,  z)/dx\  across  the 
wavefront,  which  is  physically  reasonable.  Inserting  the 
purely  spherical  form  <(>(x,  z)  =  irx‘//?(z)X  will  confirm 
that  the  integral  in  (17)  converges  to  the  expected  value 
R(z). 

To  recapitulate,  we  have  shown  that  any  real  or  non¬ 
ideal  laser  beam  with  a  wrinkled  or  nonspherical  wave- 
front  will  have  an  effective  radius  of  curvature  R{z)  given 
at  any  arbitrary  plane  z  by  the  expression  in  (9).  The  axial 
variation  of  this  effective  radius  R{z)  is  then  given  rigor¬ 
ously  by  the  Gaussian-beam-like  formula  of  (10).  Extract¬ 
ing  this  curvature  from  the  wavefront  at  any  plane  both 
converts  the  beam  profile  to  a  waist  at  that  plane  and  min¬ 
imizes  the  far-field  spread  of  the  beam  starting  from  that 
plane.  The  close  identification  between  the  propagation 
formulas  for  the  real-beam  values  of  a^(z)  and  R(z)  given 
in  (7)  and  (10)  and  the  corresponding  formulas  for  an  ideal 
TEMoo  Gaussian  beam  is  both  remarkable  and  useful . 

As  a  final  note,  we  can  observe  from  (15)  above  that 


once  the  average  phase  curvature  R(z)  has  been  extracted 
from  the  wavefront  at  any  plane  z,  the  resulting  value  of 
a]  after  the  lens  becomes  the  minimum  value  for  that  plane 
given  by 

<^s.min(2)  = 

But  this  says  that  the  total  far-field  divergence  of  the  beam 
before  the  curvature  is  extracted  can  be  written  as 

+  (^)  •  (19) 

In  other  words,  the  far-field  angular  beam  spread  a,  of  an 
arbitrary  beam  at  any  plane  z  can  be  viewed  as  the  rms 
combination  of  a  correctable  or  removable  portion 
a(z)/\R(z),  which  is  proportional  to  the  average  wave- 
front  curvature  l//?(z),  plus  an  irreducible  portion 
which  cannot  be  removed  with  any  simple  fo¬ 
cusing  element.  A  direct  analog  to  this  far-field  expres¬ 
sion  which  applies  instead  to  the  time  duration  and  the 
spectral  bandwidth  of  chirped  pulse  signals  in  time  was 
pointed  out  some  time  ago  by  Pohlig  [6] . 
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Choice  of  Clip  Levels  for  Beam  Width  Measurements 

Using  Knife-Edge  Techniques 

Anthony  E.  Siegman,  Fellow,  IEEE,  M.  W.  Sasnett,  and  T.  F.  Johnston,  Jr. 


AAsfnici— The  scanning  knife-edge  technique  provides  one  of 
the  better  methods  for  measuring  the  width  of  an  optical  beam 
having  an  arbitrary  or  irreguiar  transverse  profile.  To  imple¬ 
ment  thb  method,  however,  one  must  select  both  a  clip  level 
and  a  scale  factor  for  converting  the  measured  clip  width  into 
an  effective  width  of  the  laser  beam.  We  show  that  the  pre¬ 
ferred  clip  level  for  beam  width  measurements  using  this  tech¬ 
nique  should  be  chosen  in  the  range  between  t  »  8.5%  and  € 
«  11.6%.  With  a  suitably  matched  clip  level  and  scale  factor 
within  this  range,  the  conversion  from  measured  clip  width  to 
standard  deviation  can  be  made  exactly  accurate  for  TEMm 
Gaussian  beams,  and  the  conversion  factor  will  become  only 
slightly  inaccurate  for  a  wide  range  of  other  higher-order  or 
poorer-quality  beam  profiles. 

I.  Introduction 

Determining  the  spot  size  of  a  laser  beam,  whether 
in  the  near-  or  far-field  region  is  a  fundamental  prob¬ 
lem  in  laser  diagnostics.  The  transverse  profiles  of  real 
laser  beams  are  sometimes  veiy  irregular  in  shape,  dis¬ 
playing  multiple  peaks  and  lacking  clearly  defined  edges. 
Even  defining  the  width  of  a  laser  beam  in  a  meaningful 
way  can  be  difficult  for  an  aberrated  or  multimode  laser 
beam  with  an  irregular  intensity  profile,  let  alone  measur¬ 
ing  this  width  experimentally. 

One  precisely  defined  width  parameter  that  is  of  fun¬ 
damental  signihcance  for  any  laser  beam  is  the  standard 
deviation  evaluated  in  the  x  transverse  direction  over 
the  beam  intensity  profile.  The  variances  a]  and  al  eval¬ 
uated  in  two  transverse  directions  across  an  optical  beam 
profile  play  a  particularly  significant  role  in  theoretical 
beam  propagation  analysis  and  in  definitions  of  laser  beam 
quality.  Making  a  direct  measurement  of  the  standard  de¬ 
viations  a.,  and  for  a  real  beam  profile  can  be  difficult 
in  practice,  however,  since  this  requires  first  an  accurate 
measurement  of  the  two-dimensional  intensity  profile  / (x, 
y)  and  then  accurate  calculation  of  the  first  and  second 
moments  of  this  profile. 

A  much  simpler  and  more  widely  used  measurement 
technique  is  the  scanning  knife-edge  method.  In  this 
method,  one  scans  an  opaque  knife  edge  transversely 
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across  the  beam  profile  in,  say,  the  x  transverse  direction 
while  detecting  the  transmitted  intensity  reaching  a  large- 
area  photodetector  as  the  knife  edge  opens  up  the  beam. 
The  first  derivative  of  this  detected  intensity  versus  scan 
position  X  gives  an  indication  of  the  beam  intensity  profile 
in  the  x  direction,  integrated  over  the  y  direction. 

Although  the  complete  knife-edge  profile  contains  use¬ 
ful  information,  in  many  cases  one  wants  only  a  single 
number  to  characterize  the  effective  size  of  the  beam  in 
each  of  the  x  and  y  directions.  In  this  case,  one  can  use 
the  “clip  width”  of  the  beam  measured  in  each  of  the  two 
transverse  dimensions.  The  clip  width  for  the  beam  in 
either  direction  is  defined  as  the  width  between  points  in 
the  scan  where  the  transmitted  intensity  passing  the  knife 
edge  rises  from  a  lower  clip  level  ePo  to  an  upper  clip 
level  (1  -  t)Po  where  Pq  is  the  total  power  in  the  beam. 
Clip  levels  in  the  range  around  10%  and  90%  are  often 
used. 

Selecting  an  optimum  clip  width  and  an  associated  scale 
factor  for  such  knife-edge  measurements  is  the  primary 
topic  of  this  paper.  We  show  that  1)  there  is  an  unavoid¬ 
able  error  or  uncertainty  in  converting  such  clip  width 
measurements  into  more  theoretically  significant  mea¬ 
sures  such  as  the  standard  deviation  of  the  beam  in  the 
same  direction;  and  2)  any  clip  level  e  in  the  range  be¬ 
tween  about  8.5  %  and  11.6%,  together  with  a  scale  factor 
based  on  an  ideal  TEMqo  Gaussian  beam,  can  be  about 
equally  valid  as  a  standard  definition  for  knife-edge  mea¬ 
surements  on  a  wide  variety  of  real  laser  beams. 

II.  Clip  Level  and  Scale  Factor  Selection 


The  clip  width  for  any  specific  beam,  measured 
using  a  specified  clip  level  e,  is  of  course  a  definite  and 
fairly  easily  measured  quantity.  The  standard  deviation 
of  the  beam  in  the  knife-edge  scanning  direction,  which 
is  a  quantity  of  more  theoretical  significance,  can  then  be 
estimated  from  the  measured  clip  width  using  a  suitable 
scaling  factor  S  such  that  a,  =  D^/S.  There  is,  however, 
an  unavoidable  uncertainty  in  going  from  measured  clip 
width  to  the  more  fundamental  standard  deviation  in  this 
fashion.  The  exact  scale  factor  or  the  exact  value  of  the 
ratio  Dr/ a,  for  an  arbitrary  beam  will  depend  not  only  on 
the  clip  level  employed  but  also  to  a  certain  extent  on  the 
exact  shape  of  the  laser  beam. 


A.  Scale  Factor  Errors  in  Clip  Level  Measurements 
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To  put  this  in  another  way,  different  beams  having  the 
same  standard  deviation  a,  will  give  slightly  different 
measured  clip  widths  depending  on  the  exact  profile 
and  modal  content  of  the  beam,  which  is  assumed  to  be 
unknown.  Hence,  there  will  always  be  some  uncertainty 
in  the  scale  factor  S  that  should  be  used  to  transform  from 
the  measured  clip  width  D^.  to  the  more  theoretically  in¬ 
teresting  value  (T,.  In  practice,  one  selects  the  scaling  fac¬ 
tor  S  to  match  the  calculated  value  D,./a,  for  some  stan¬ 
dard  beam  profile,  preferably  a  TEMoo  Gaussian  profile. 
This  scale  factor  will  then,  in  general,  exhibit  some  error 
when  applied  to  other  beam  profiles.  Our  primary  objec¬ 
tive  in  this  paper  is  to  show  how  to  estimate  and  minimize 
these  errors  for  a  wide  range  of  real  beam  profiles,  by 
examining  the  theoretical  clip  level  versus  clip  width 
curves  for  a  representative  sample  of  different  beam  pro¬ 
files. 

B.  Low-Order  Gaussian  Beams 

Fig.  1  shows  by  way  of  example  the  fractional  power 
falling  outside  a  given  clip  width  on  one  side  of  the 
beam  only  (which  is  the  same  thing  as  the  clip  level  f )  for 
several  different  Hermite-Gaussian  (HG)  or  Lagueire- 
Gaussian  (LG)  optical  beam  profiles.  The  fractional  power 
in  each  case  is  plotted  versus  the  clip  width  normalized 
to  the  standard  deviation  for  that  specific  beam.  The 
four  mode  profiles  considered  in  this  figure  are;  1)  the 
lowest-order  or  n  =  0  Hermite-Gaussian  mode,  which 
corresponds  also  to  a  TEMqo  mode  in  two  transverse  di¬ 
mensions;  2)  the  first  antisymmetric  orn  =  I  HG  mode, 
which  would  correspond  to  a  TEM|o  but  not  an  TEMqi 
HG  mode  in  two  transverse  dimensions;  3)  the  “donut” 
mode,  or  Laguerre-Gaussian  TEMqi  mode,  which  often 
appears  as  a  perturbation  in  stable-resonator  laser  oscil¬ 
lators,  and  which  corresponds  also  to  an  incoherent  su¬ 
perposition  of  TEM|o  and  TEMqi  HG  modes;  and  finally 
4)  an  annular  ring  mode  which  represents  an  asymptotic 
approximation  for  a  very  high-order  HG  mode,  as  dis¬ 
cussed  below.  Note  again  that  the  horizontal  axis  in  each 
case  is  scaled  to  the  standard  deviation  of  each  given 
mode,  and  not  to  a  common  Gaussian  spot  size  w  char¬ 
acterizing  the  family  of  modes. 

The  dashed  lines  on  this  plot  show  that  if  one  selects, 
for  example,  a  clip  level  e  =  11.59%  and  a  scale  factor 
5  =  2.392,  knife-edge  measurements  using  these  param¬ 
eters  will  yield  exactly  accurate  values  of  cr,  for  beams 
consisting  of  either  pure  TEMqq  or  pure  HG  n  -  I  (al¬ 
though  not,  in  general,  for  mixtures  of  these  two  modes). 
Using  this  same  clip  level  and  scale  factor  for  the  other 
modes  shown  will  slightly  underestimate  the  scale  factor 
or  overestimate  the  value  of  a,,  with  an  error  in  the  range 
of  5-10%.  because  the  scale  factor  S  will  be  slightly  dif¬ 
ferent  than  the  exact  Dja^  values  for  these  modes  at  this 
selected  clip  level.  (Note  again  that  in  a  simple  automated 
knife-edge  measurement,  only  the  clip  level  e  and  clip 
wid'h  will  be  directly  determined,  and  one  will  in  gen¬ 
eral  ha^e  no  direct  knowledge  of  the  exact  beam  profile 
being  measured.) 


Fig.  I.  Fractional  power  falling  outside  a  given  clip  width  D,.  on  one  side 
plotted  versus  O,  /o,  for  several  low-order  Hermite-Gaussian  modes.  Shown 
in  this  figure  are  the  two  lower-order  Hermite-Gaussian  modes  n  =  0  and 
n  =  1;  the  donut  mode  (incoherent  mixture  of  TEM,n  and  TEMoi);  a  ring 
mode  which  corresponds  to  the  asymptotic  Hermite-Gaussian  case  with  n 
-•  oo;  and  two  suggested  clip  levels  at  r  =  8.5%  and  t  =  11.6%.  The 
insert  shows  experimental  data  for  the  TEM,x>  and  donut  modes  measured 
on  the  same  argon-ion  laser  using  a  knife-edge  system  with  an  adjustable 
clip  level. 

From  this  figure  it  appears  that  a  useiul  alternative 
choice  could  be  to  select  a  clip  level  e  =  8.485%  and  a 
scale  factor  S  =  2.746.  This  choice  would  give  exact  re¬ 
sults  for  both  the  TEMqo  mode  and  the  donut  mode,  with 
scale-factor  errors  of  plus  or  minus  a  few  percent  for  the 
other  modes  shown.  As  an  experimental  confirmation  of 
this  point,  the  smaller  plot  in  Fig.  1  displays  clip  widths 
measured  on  the  same  argon-ion  laser  beam  under  two 
different  operating  conditions,  using  a  scanning  knife-edge 
system  with  an  electronically  adjustable  clip  level.  In  one 
case,  the  laser  was  operated  in  a  TEMqq  mode  as  deter¬ 
mined  by  a  suitable  intracavity  mode-control  aperture.  In 
the  other  case,  a  Brewster  plate  with  a  small  defect  burned 
into  one  spot  on  the  plate  was  positioned  within  the  ion 
laser  cavity  so  as  to  force  oscillation  in  the  donut  mode 
within  the  same  cavity  mirrors.  Single-transverse-mode 
operation  was  confirmed  in  each  case  by  the  absence  of 
any  transverse  mode  beats.  The  measured  clip  widths  in 
the  two  cases  were  then  scaled  to  the  appropriate  standard 
deviations  assuming  a  common  Gaussian  spot  size  pa¬ 
rameter  w  for  the  two  modes  in  the  same  resonator,  with 
w  being  measured  on  the  TEMqq  mode.  The  measured  clip 
level  curves  intersect  at  almost  exactly  8.5%  as  predicted, 
and  this  intersection  occurs  at  a  D^fa^  ratio  of  2.75  as 
also  predicted  by  the  theory. 

Either  of  the  two  clip  levels  e  =  8.5%  or  f  =  1 1.65% 
might  be  regarded  as  “natural”  choices,  at  least  for  beams 
made  up  of  HG  mode  mixtures.  We  will  show  more  gen¬ 
erally,  however,  that  any  choice  of  clip  level  between 
these  two  extremes  will  give  a  very  similar  range  of  errors 
for  a  very  wide  range  of  optical  beam  profiles,  whether 
Hermite-Gaussian  in  character  or  not. 

C.  Higher-Order  Hermite-Gaussian  Beams 

Fig.  2  shows  curves  of  the  same  structure  as  Fig.  1  with 
the  donut  mode  removed  for  simplicity  and  with  higher- 
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Fig.  2.  Clip  level  plot  siinilar  to  Fig.  I  with  the  donut  mode  removed  and 
higher-order  HG  modes  from  n  =  2  to  n  =  7  added. 

order  Hermite-Gaussian  modes  with  indexes  from  n  =  2 
up  to  n  =  7  added.  From  this  figure  it  is  evident  from  the 
curves  for  the  higher-order  HG  modes  oscillate  in  a  reg¬ 
ular  fashion  about  the  “ring  mode”  or  asymptotic  HG  n 
00  curve  introduced  in  Fig.  1.  This  behavior  can  be 
understood  as  follows. 

The  one-dimensional  intensity  profiles  for  higher-order 
Hermite-Gaussian  beams  have  the  general  form  l„(x)  = 
Hl(y/2x  /w)  exp  (-2jt^/tv*),  with  standard  deviation 

=  s/2n  +  1  w/2.  The  Hermite-Gaussian  intensity  profile 
for  the  specific  case  n  ~  10  is  shown  by  the  solid  line  in 
Fig.  3.  But  this  mode  profile  is  also  the  well-known  quan¬ 
tum  probability  distribution  for  a  one-dimensional  quan¬ 
tum  harmonic  oscillator.  As  many  quantum  theory  texts 
point  out,  when  the  mode  index  n  becomes  large,  these 
Hermite-Gaussian  profiles  with  the  high-frequency  peri¬ 
odic  variation  averaged  out  approach  the  one-dimensional 
probability  di.stribution  for  a  classical  particle  executing 
sinusoidal  harmonic  motion  in  one  dimension.  The  latter 
distribution  has  the  normalized  form 


(1) 

with  the  outer  width  and  the  standard  deviation  being  re¬ 
lated  by  JCoo  =  7/20,.  This  limiting  profile  is  indicated  by 
the  dashed  line  in  Fig.  4. 

But  a  classical  particle  oscillating  sinusoidally  in  one 
transverse  dimension  has  the  same  probability  distribution 
with  respect  to  that  dimension  as  does  a  particle  traveling 
at  constant  velocity  around  a  circular  orbit  in  two  trans¬ 
verse  dimensions.  As  a  result,  the  asymptotic  one-dimen¬ 
sional  intensity  distribution  for  HG  modes  as  n  -»  00  also 
represents  equally  well  the  one-dimensional  distribution 
for  an  infinitely  thin  annular  “ring”  beam  with  equal  in¬ 
tensity  around  the  annulus. 

D.  Other  Elementary  Beam  Profiles 
The  focus  on  Hermite-Gaussian  modes  in  Figs.  I  and 
2  may  merit  further  discussion.  The  mode  patterns  for  real 
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Fig.  3.  Intensity  profile  l„ix)  for  the  Hermite-Gaussian  mode  with  n  =  10 
(solid  line),  and  the  asymptotic  form  of  this  profile  for  n  -•  oo  (dashed 
line). 


Fig.  4.  Clip  level  plot  similar  to  Fig.  1  with  various  other  simple  beam 
profiles  added. 

laser  beams  will  in  many  cases  not  be  single  Hermite- 
Gaussian  or  Laguerre-Gaussiun  modes,  but  rather  may 
have  a  variety  of  more  complex  forms.  These  more  com¬ 
plex  profiles  can  nonetheless  be  viewed,  if  one  wishes,  as 
either  coherent  or  incoherent  superpositions  of  Hermite- 
Gaussian  modes  (other  basis  sets  could  of  course  be  used 
equally  as  well).  Our  assumption  is  that  if  one  uses  a  mea¬ 
surement  procedure  that  is  demonstrably  accurate  for  the 
limiting  case  of  an  n  =  0  Hermite-Gaussian  beam,  and 
that  has  a  small  nnd  bounded  error  for  all  higher-order 
Hermite-Gaussian  modes,  then  this  procedure  1 )  will  be 
most  accurate  for  the  most  important  case,  namely,  low- 
urder  and  thus  high-quality  beam  profiles;  and  2)  it  may 
have  only  small  and  well-understood  inaccuracies  for  ar¬ 
bitrary  beams,  which  can  always  be  described  as  mixtures 
of  higher-order  Hermite-Gaussian  modes.  The  latter  as¬ 
sumption  is.  unfortunately,  not  entirely  correct,  as  we  will 
demonstrate  in  the  following. 

To  supplement  the  HG  and  LG  examples  of  Figs.  1  and 
2,  several  other  elementary  beam  profiles  including  uni¬ 
form  top  hat  and  unifvJim  rectangular  beams,  and  dia¬ 
mond-shaped  and  “picture-frame”  beam  profiles,  have 
been  plotted  in  Fig.  4  along  with  the  same  low-order 
Gaussian  beams  as  in  Fig.  1.  The  picture-frame  mode, 
representing  a  beam  profile  in  the  form  of  a  very  thin 
square  annulus,  is  something  of  an  “odd  man  out”  in  this 
plot,  with  a  clip  level  that  jumps  instantaneously  from  0 
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to  25%  at  the  point  where  the  knife  edge  encounters  the 
flat  edge  of  the  picture  frame. 

Fig.  5  shows  an  expanded  view  of  the  clip  level  region 
between  6  =  7%  and  £  =  13%  with  all  the  curves  shown 
in  Figs  1-4  included.  Three  potential  clip  levels  and  scale 
factors  are  also  indicated,  namely,  the  clip  level  at  11.6% 
where  the  TEMqo  and  HG  «  =  1  modes  will  have  equal 
Dj  / Oj  values;  the  clip  level  at  8.5%  where  the  TEMqo  and 
the  donut  modes  will  have  equal  £>,  /a,  values;  and  a  com¬ 
promise  value  arbitrarily  located  at  the  10%  clip  level. 
The  sum  total  of  the  results  in  Figs.  1-5  give  a  reasonable 
(although  unfortunately  not  complete)  indication  of  the 
spread  in  scale  factors  that  might  be  encountered  for  real 
laser  beams  in  practice.  The  middle  choice  e  =  10%  might 
represent  a  reasonable  compromise  between  the  two  outer 
values,  with  a  range  of  scale  factors  which  remains  small 
for  low-order  Hermite-Gaussians  and  is  distributed  more 
or  less  symmetrically  about  the  TEMqq  scale  factor. 

In  all  of  these  cases,  we  chose  a  scale  factor  S  which 
will  be  exactly  correct  for  a  TEMqq  Gaussian  beam.  This 
is  done  for  two  reasons.  First,  in  real  lasers  using  stable 
resonators  and  having  good  beam  quality,  the  output  mode 
will  often  be  very  close  to  a  TEMqq  beam;  to  get  the  most 
accurate  results  in  this  case,  one  should  select  a  scale  fac¬ 
tor  which  applies  to  that  beam.  Second,  in  the  more  gen¬ 
eral  case  of  more  irregular  beam  profiles,  the  "beam  qual¬ 
ity”  for  an  arbitrary  beam  is  most  often  defined  in  a 
manner  which  uses  the  TEMqq  beam  as  a  reference  stan¬ 
dard.  Beam  quality  measurements  using  the  TEMqq  scale 
factor  will  therefore  become  most  accurate  in  the  impor¬ 
tant  limit  of  good  quality  beams,  and  the  unavoidable  un¬ 
certainty  in  the  scale  factor  will  only  show  up  in  the  less 
important  case  of  poorer  quality  beams. 

E.  Pedestal  Beams 

Despite  the  reasonable  range  of  error  values  indicated 
by  Figs.  1-5,  one  should  be  aware  that  there  can  always 
be  extreme  beam  profiles  having  very  different  values  of 
both  higher  and  lower,  than  the  ranges  indicated 
for  a  given  choice  of  clip  level  in  Figs.  1-5.  Consider, 
for  example,  a  beam  profile  consisting  of  two  outer  spikes 
or  delta  functions  located  at  jt  =  ±  Ma  with  A/  >  1 ,  each 
containing  just  under  e  of  the  total  beam  power,  and  two 
spikes  at  jr  =  ±a  each  containing  0.5  -  £  of  the  beam 
power.  The  normalized  clip  width  for  chip  level  £  in  this 
case  will  be  D^./a^  —  2/yJ\  +  2(M~  —  1)£.  which  can 
be  made  as  small  as  desired  for  M  »  \.  Alternatively, 
consider  a  beam  profile  with  two  outer  spikes  at  .x  =  ±a 
containing  just  over  £  of  the  total  power  each,  with  all  the 
remaining  energy  in  a  narrow  spike  at  jr  =  0.  The  nor¬ 
malized  clip  width  in  this  case  will  be  given  by  D,  fa,  = 
l/£,  which  will  be  much  larger  than  any  of  the  values 
shown  in  Figs.  1-5. 

Another  bothersome  case  is  the  simple  example  of  the 
far-field  pattern  of  a  single  slit,  which  has  the  general  form 
(at  lea,st  in  the  elementary  paraxial  approximation)  of  l(x) 
=  sin'  x/.x^.  The  problem  in  this  case  is  that,  because  of 
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Fig.  5.  Expanded  view.';  of  selected  clip  level  curves  from  Figs.  1-5  over 
the  range  f  =  7%toe=  13%.  Three  suggested  clip  levels  and  the  asso¬ 
ciated  scale  factors  for  the  TEMoo  mode  are  indicated. 

the  very  slow  1/x'  falloff  of  this  profile  at  large  x,  the 
value  of  <7,  actually  diverges  logarithmically.  Hence,  the 
value  of  D,  /a,  cannot  even  be  defined  without  putting  in 
some  truncation  of  the  sin^  x/x^  function  at  some  large 
but  finite  value  of  X. 

Still  another  not  uncommon  example  would  be  a  "ped¬ 
estal  beam,"  which  might  be  modeled  as  consisting  of  a 
TEMqq  Gaussian  primary  beam  with  a  Gaussian  spot  size 
w,  containing  a  fraction  1  -  a  of  the  total  power,  super¬ 
imposed  on  a  "pedestal"  modeled  as  a  second  TEMqq 
beam  with  a  spot  size  hs  containing  the  remaining  a  of 
the  total  beam  power.  High-power  laser  oscillators  will 
sometimes  produce  beams  of  this  general  character,  with 
a  central  beam  lobe  of  rather  good  quality  containing  some 
fraction  of  the  total  beam  pov^  ’r,  surrounded  by  a  broad 
low  pedestal  containing  the  remaining  power.  (Such 
beams  are  in  fact  of  quite  poor  beam  quality,  despite  their 
well-formed  central  lube.) 

Figs.  6  and  7  sliuw  clip  level  plots  for  a  series  of  such 
pedestal  beams;  first  for  a  fixed  pedestal  width  ratio  M  = 
us/w’i  and  varying  fractional  power  a,  and  then  for  fixed 
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Fig.  6.  Clip  level  plot  similar  to  Figs.  1.2.  and  4  fora  ' 'pedestal’'  beam 
in  which  the  pedestal  Gaussian  U  M  =  3  times  wider  than  the  central 
Gaussian,  with  varying  fractions  a  =  20.  40.  60.  and  80%  of  the  total 
beam  power  in  the  pedestal. 


Fig.  7.  Clip  level  curves  for  pedestal  beams  with  5%  of  the  total  beam 
power  in  the  pedestal,  with  the  pedestal  Gaussian  being  a  varying  factor  M 
=  2.  3.  4.  and  5  times  wider  than  the  central  Gaussian. 


pedestal  power  a  and  varying  width  ratio  M.  The  general 
conclusion  is  that  over  most  of  the  clip  level  range,  these 
beams  have  a  D^/ a,  ratio  that  falls  below  the  TEMqo 
value;  and  the  scale  factor  errors  in  measuring  these  beams 
can  become  quite  sizable  for  wide  pedestals  or  pedestals 
that  contain  any  significant  fraction  of  the  total  power. 
There  is  probably  no  simple  way  around  this  difficulty, 
except  the  practical  answer  of  using  a  spatial  filter  in  the 
real  laser  beam  to  eliminate  most  of  the  pedesta',  and  then 
measuring  the  remaining  output  beam  after  the  spatial  fil¬ 
ter.  In  any  event,  the  scale  factor  error  for  a  pedestal  beam 
becomes  sizable  only  for  rather  poor  quality  beams. 

F.  Mixed-Mode  Gaussian  Beams 

Finally,  one  might  guess  that  if  a  clip  level  is  chosen 
at  which  two  different  HG  modes  have  the  same  scale  fac¬ 
tors  (as  at  6  =  1 1.6  or  8.5%),  then  an  incoherent  mixture 
of  these  two  modes  might  have  the  same  D,  /a,  value  also. 
This  turns  out  to  be  nearly  but  not  exactly  true.  If  one 


(b) 

Fig.  8.  Fractional  changes  in  the  clip  width  to  standard  deviation  ratio 
D,  /o,  versus  fractional  mode  mixture  for  incoherent  superpositions  of  La- 
guerre-Gaussian  modes.  Results  are  shown  for  clip  levels  (a)  t  =  8.5% 
(where  the  (X)  and  01*  modes  have  the  same  scale  factor)  and  (b)  t  = 
10.0%.  In  each  segment  of  the  plots,  the  beam  profile  is  an  incoherent 
mixture  of  the  two  adjacent  Laguerre-Gaussian  mode  profiles,  with  the 
fractional  percentage  varying  linearly  along  the  horizontal  axis  as  one  moves 
from  one  mode  pattern  to  the  next. 


considers  a  beam  profile  consisting  of  an  incoherent  mix¬ 
ture  of  two  modes  with  varying  fractions  of  each  mode, 
the  clip  width  and  the  standard  deviation  for  the  mixed 
mode  will  both  change  in  a  slightly  different  manner  as 
the  mixture  ratio  is  changed;  and  the  ratio  will  also 

change  slightly  even  if  it  returns  to  the  same  value  at  both 
end  points.  To  illustrate  this  point.  Fig.  8(a)  and  (b)  shows 
the  fractional  change  in  the  0,  /a,  ratio,  relative  to  the 
TEMqo  Gaussian  mode,  for  various  incoherent  binary 
mixtures  of  modes  at  the  two  different  clip  levels  e  = 
1 1 .6%  and  €  =  10%.  In  each  segment  of  these  plots,  the 
beam  profile  is  assumed  to  be  an  incoherent  mixture  of 
the  two  adjacent  HG  mode  profiles,  with  the  fractional 
percentage  of  each  mode  varying  linearly  as  one  moves 
from  one  intersection  point  to  the  next.  The  general  con¬ 
clusion  is  that  there  is  a  small  variation  of  the  D,  /a,  ratio 
with  mode  mixture  ratio,  even  when  the  endpoints  have 
the  same  limiting  value.  This  variation  is  in  general  so 
small,  however,  that  it  will  have  negligible  importance  in 
affecting  the  accuracy  of  any  practical  measurement. 

111.  Summary 

The  curves  in  Figs.  1-8  should  give  a  good  indication 
of  the  variation  in  Djo^  likely  to  be  encountered  with 
real  laser  beams,  and  thus  provide  guidance  in  selecting 
a  clip  level  for  standardized  knife-edge  measurements. 
The  primary  conclusion  of  this  paper  is  that  any  cUp  level 
in  the  range  between  f  =  8.5%  and  f  =  1 1 .6%  could  be 
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about  equally  useful  as  a  standard  value  for  knife-edge 
beam-width  measurements.  The  exact  clip  level  can  be 
chosen  based  on  other  considerations  such  as  numerical 
simplicity  or  experimental  convenience.  Smaller  clip  lev¬ 
els  are  probably  less  convenient  experimentally,  from 
considerations  of  noise  and  detector  sensitivity.  Hence, 
either  e  =  10%  with  a  scale  factor  =  2.563,  ort  =  11.6% 
with  a  scale  factor  D^/o^  =  2.39  might  be  equally  useful. 

The  scale  factor  Dja^  to  be  used  for  converting  mea¬ 
sured  clip  widths  into  (approximate)  values  of  a,  for 
real  laser  beams  should,  in  any  case,  be  the  value  appro 
priate  a  TEMqo  Gaussian  beam.  Proposals  have  been  made 
by  some  standards  groups  to  define  the  half-width  in 
the  TT-direction  for  a  real  laser  beam  to  be  =  la^.  This 
half-width  quantity  for  real  laser  beams  would  then  match 
up  directly  with  the  1 J  e~  half-width  parameter  or  Gauss¬ 
ian  spot  size  parameter  w,  =  2a,  that  is  widely  used  for 
TEMqo  Gaussian  beams.  If  this  proposal  were  adopted, 
and  if  the  10-90%  choice  were  also  adopted  as  a  standard 
for  the  clip  level  e,  then  the  measured  full  width  D  =  2W 
for  a  laser  beam  would  be  related  to  the  measured  clip 
width  by 

We  note  again  that  in  making  a  measurement,  one  will 
usually  not  know  the  modal  content  of  the  beam,  so  that 
one  can  only  measure  the  clip  width  at  a  predetermined 
clip  level  and  then  apply  the  predetermined  scale  factor 
to  obtain  an  estimated  value  of  the  standard  deviation  or 
other  selected  beam  size  parameter.  Since  at  least  some 
real  lasers  approach  pure  TEMqq  behavior,  often  with 
small  amounts  of  TEMqi  or  TEMio  mixed  in,  this  ap¬ 
proach  will  become  exactly  correct  in  the  TEMqq  limit. 
Measurements  using  this  approach  will  then  become  grad¬ 
ually  less  accurate  for  higher-order  or  poorer-quality 
beams,  depending  upon  the  exact  beam  profile. 

There  will  always  be  extreme  beam  profiles  for  which 
the  error  in  converting  from  clip  width  to  standard  devia¬ 
tion  can  become  rather  large.  There  appears  to  be  no  easy 
way  around  this  problem.  One  can,  of  course,  gain  a  sub¬ 
stantial  amount  of  additional  information  by  using  the 
knife-edge  technique  to  record  the  complete  profile  of 
transmitted  intensity  versus  position,  and  then  using  this 
information  to  determine  the  character  of  the  mode  profile 
more  accurately.  But  this  requires  considerably  more  data 
handling  and  recording  than  the  simple  procedure  of  re¬ 
cording  a  single  clip  width  D,.  at  a  single  standardized  clip 
level  f. 


Appendix  A 

Clip  Level  Forumulas 

The  analytical  clip  level  formulas  for  the  modes  dis¬ 
cussed  in  this  paper  are  summarized  in  this  Appendix. 

Given  an  arbitrary  transverse  intensity  profile  Hx.  y), 
the  total  beam  profile  P,,  and  the  standard  deviation  o,  in 


the  jr-direction  for  an  optical  beam  are  given  by 

SOo  ^00 

j  I{x,  v)  dx  dy  and 

—  00  Joo 

SQo  A  oo 

1  x-l(x,  y)  dx  dy.  (2) 

—  oo  J  00 

The  clip  width  in  the  jr-direction  for  a  given  clip  level  e 
will  then  be  given  by  D,.  =  —  O),  where  a,  and  oj  are 

the  points  at  which 


dy  /(jt.  y)  =  e  Pq  and 

dy  /(.c,  y)  =  (1  -  e)Po. 


(3) 


If  we  consider  only  symmetric  beam  profiles  for  simplic¬ 
ity,  then  fli  =  02  =  A  /2  and  the  clip  level  e  can  be  re¬ 
lated  to  the  clip  width  D,.  by 

I 

e  =  3 


•  lyi  /  -  » 

-  dx 

J-D,/:  Jo 


dy  I(x,  y) 


(4) 


The  knife-edge  method  of  course  automatically  integrates 
over  the  y  variation  when  the  knife  moves  in  the  T:-direc- 
tion,  and  we  can  drop  the  y  variation  for  any  cases  where 
/(jc,  y)  is  separable  in  the  form  I{x,  y)  =  /,(jt)  x  /,(  y). 
Clip  level  formulas  for  the  beams  considered  in  this  paper 
are  then  given  by  the  following  formulas. 

I)  Lowest-Order  Hermite-Gaussian  Mode:  The  clip 
level  for  the  lowest-order  HGq  mode  profile  l^ix)  =  exp 
(-2jc'/h'’),  or  for  the  corresponding  TEMqq  or  TEMqq 
HG  modes,  is  given  by 
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\ 
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(5) 


2)  Next  Higher-Order  Hermite-Gaussian  Mode:  The 
clip  level  for  the  HG  n  =  1  mode  or  for  the  corresponding 
TEM|o,  but  not  TEMo,,  mode  is  given  by 
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(6) 


3)  The  "Donut"  Mode:  This  mode,  which  is  often 
seen  in  stable-resonator  lasers  in  which  the  mode-con¬ 
trolling  aperture  is  made  slightly  too  large,  can  be  viewed 
either  as  a  Laguerre-Gaussian  TEM(*  mode  with  profile 
given  by  l{r.  6)  =  r'  exp  l-2rV‘'''l.  ot  as  the  inco¬ 
herent  superposition  of  equal  amounts  of  the  TEM,o  and 
TEMq,  HG  modes  with  profiles  given  by  /|o  =  x'  exp 
\-2(x"  +  y')/M  "l  and /|,|  =  y"  exp  ( -2(.v-  -(-  y’)/H  ’|. 
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In  either  case,  the  clip  level  is  given  by 


e 


2 


+ 


1 

2 a, 


exp 


(7) 

4)  ‘  ‘Ring  ’  ’  Mode  {Asymptotic  Hermite-Gaussian 
Mode):  As  discussed  in  the  text,  the  asymptotic  limit  for 
an  HG  mode  as  n  -*  oo  has  the  same  properties  in  one 
transverse  dimension  as  does  a  narrow  ring  annulus  beam 
in  two  dimensions.  The  D/a,  ratio  as  a  function  of  clip 
level  for  either  case  is  given  by 


e 


(8) 


5)  “Top  hat”  Mode:  A  circular  uniform-amplitude  top 
hat  beam  has  clip  level  and  clip  width  related  by 


t 


2 

2 


(6  A  sin  6  cos  6) 


where  0  =  sin  '  (D/2a,). 


(9) 


6)  Uniform  Slit  or  Rectangular  Mode:  A  uniformly  il¬ 
luminated  slit,  or  any  form  of  equivalent  rectangular  beam 
scanned  across  either  of  its  principal  axes,  will  have  the 
form 


6 


(10) 


7)  “Diamond”  Mode:  A  "diamond”  profile  (a  ui’i- 
formly  illustrated  square  beam  rotated  by  45°)  obeys  the 
formula 


€ 


(II) 


8)  '  ‘Picture  Frame ' '  Mode:  A  very  thin  square  annu¬ 
lus  (the  square  analog  of  the  ring  mode)  obeys  the  formula 

D  (Ve  X  (2  ~  4(),  0.25  <  €  <  0.5 

-  =  r-  (12) 

C’^6.  e  <  0.25. 

9)  Higher-Order  Hermite-Gaussian  Modes:  Analytic 
formulas  can  be  developed  if  neces.sary  for  higher-order 
HG  modes  with  n  >  2.  but  in  most  cases  it  proves  simpler 


to  evaluate  the  properties  of  these  modes  by  numerical 
evaluation  using  a  computer  package  such  as  Matljcma- 
tica,  and  this  was  done  to  calculate  the  results  in  Figs.  2 
and  5. 

10)  Pedestal  Beams:  Consider  a  beam  consisting  of  a 
central  TEMoo  Gaussian  with  spot  size  w,  containing  1  - 
a  of  the  total  power,  and  a  pedestal  with  spot  size 
containing  the  remaining  fraction  a  of  the  total  power. 
The  clip  level  is  then  given  by 
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Single  Pulse  Laser  Beam  Quality  Measurements  Using  a  CCD  Camera  System 
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We  have  assembled  a  CCD-camera-based  system  that  is  capable  of  measuring  laser  beam 
quality  on  a  single-shot  pulsed  basis.  The  system  has  been  tested  using  a  pulsed  HeNe  laser 
beam  with  added  spatial  distortion. 

The  concept  of  beam  quality  as  a  measure  of  the  multi-mode  or  "times-diffraction- 
limited"  character  of  a  laser  beam  has  been  of  interest  since  the  early  days  of  the  laser  field.  ^ 

Only  recently,  however,  has  there  begun  to  be  convergence  on  standard  methods  for  defining  and 
measuring  laser  beam  quality.^  In  the  past  few  years  commercial  instruments  for  measuring 
laser  beam  quality  on  a  CW  basis  have  become  available  and  have  been  found  to  yield  useful 
new  information  on  laser  physics.^  Our  device  enables  similar  investigations  of  pulsed  lasers. 

We  have  developed  a  CCD-camera-based  system  that  makes  single-shot  beam  quality 
measurements  on  pulsed  laser  sources.  The  system  employs  three  CCD  cameras,  a  Macintosh 
computer  and  three  video  frame  grabber  boards,  as  shown  in  Fig.  1.  The  cameras  monitor  the 
transverse  beam  profile  as  the  beam  propagates  through  a  small  beam  waist  formed  by  a  high 
quality  lens.  The  frame  grabber  boards  convert,  or  digitize,  the  video  images  of  the  laser  beam 
into  arrays  of  numbers  that  the  computer  uses  to  calculate  beam  widths  at  each  plane.  The 
combination  of  the  beam  waist  size  and  the  far-field  divergence  of  the  laser  beam  provides  a 
measure  of  the  beam  quality. 

To  ensure  an  accurate  measurement  of  beam  quality  the  samples  of  the  transverse  beam 
profile  must  be  taken  near  the  beam  waist  and  several  Rayleigh  ranges  away  from  the  waist. 

Use  of  a  short  Rayleigh  range  enables  measurement  in  a  confined  laboratory.  In  our  experiments 
a  high  quality  10  cm  focal  length  lens  produces  a  very  small  waist  with  a  typical  Rayleigh 
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range  of  a  few  centimeters. 

The  transverse  profiles  are  measured  with  standard  Sony  SSC-D5  CCD  cameras.  The 
frame  grabber  boards  quantize  the  continuous  video  information  from  the  cameras  into  a  finite 
number  of  light  intensity  samples  (pixels  on  the  computer  screen).  They  also  quantize  the 
intensity  at  each  sample  into  a  finite  number  of  levels.  To  reduce  quantization  error,  each 
transverse  beam  image  should  fill  as  much  of  the  area  of  its  camera's  CCD  array  as  possible 
without  overflow.  Also,  the  intensity  should  be  as  great  as  possible  without  saturating  the 
camera.  If  only  10%  of  the  maximum  range  of  intensity  levels  is  used,  the  measured  beam  width 
can  be  in  error  by  as  much  as  5%.  Below  that,  the  measured  width  falls  quickly  with  intensity. 

One  CCD  camera  coupled  with  a  microscope  objective  monitors  the  transverse  profile  of 
the  laser  beam  at  the  small  waist  formed  by  the  focusing  lens  (Fig.  1).  The  microscope  objective 
allows  measurement  of  waists  as  small  as  a  few  microns.  The  distance  between  the  camera  and 
the  microscope  objective  is  fixed  so  that  objectives  can  be  easily  interchanged  and  the 
magnification  simply  read  off  the  objective  in  use.  As  mentioned  above,  the  power  of  the 
objective  must  be  chosen  to  project  the  beam  waist  onto  the  largest  possible  area  of  the  CCD 
detector.  For  example,  with  256  intensity  levels  and  a  near-gaussian  laser  beam,  the  gaussian 
spot  size  should  be  roughly  30%  of  the  short  edge  of  the  camera's  CCD  array.  At  this  size,  with 
the  intensity  at  the  center  of  the  beam  just  below  saturation,  the  intensity  at  the  edge  of  the 
camera  will  be  below  the  lowest  quantization  level. 

The  exact  position  of  the  beam  waist  is  located  by  moving  the  camera  and  microscope 
objective  as  a  unit,  as  in  a  video  microscope'^,  and  making  measurements  on  repeated  shots  of  the 
pulsed  laser,  looking  for  a  minimum  spot  size.  In  the  case  of  an  astigmatic  beam,  the  two 
transverse  axes  must  be  measured  separately.  During  positioning  of  the  microscope  assembly, 
the  computer  calculates  a  continuous  rough  estimate  of  beam  width  using  a  rapid  width 
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measurement  technique ,  such  as  a  simulated  knife  scan.  A  more  accurate  width  measurement 
algorithm  is  used  during  data  collection. 

The  two  other  cameras  monitor  the  transverse  beam  profile  at  varying  distances  beyond 
the  waist  (Fig.  1).  As  the  laser  beam  spreads,  one  of  the  cameras  is  moved  to  a  position  where 
the  laser  beam  just  fills  the  entire  CCD  array,  as  mentioned  before.  The  third  camera  is  placed 
somewhat  closer,  to  obtain  a  measurement  at  a  different  distance  from  the  beam  waist.  The 
beam  profile  on  this  camera  will  not  quite  fill  the  array. 

Any  imperfect  optical  components  in  the  optical  path  will,  of  course,  change  the  beam 
quality  we  are  attempting  to  measure.  Our  design  reduces  the  number  of  these  components  in  the 
measuring  device  to  a  10  cm  focal  length  lens,  a  microscope  objective  and  two  high-quality 
beamsplitters  (Rg.  1).  The  CCD  cameras  were  used  without  lenses.  In  general,  a  variable 
neutral  density  filter  will  also  be  required  to  provide  control  of  the  over-all  intensity.  In  our 
testing  of  this  s)«tem,  however,  we  are  able  to  control  the  average  intensity  of  the  beam  by 
changing  the  pulse  length.  Separate  neutral  density  filters  are  not  required  for  each  camera  if 
the  beamsplitters  are  selected  so  that  the  maximum  (spatial)  intensity  is  just  below  saturation 
on  each  camera.  In  our  set-up  the  first  beanasplitter  reflects  roughly  33%  and  the  second 
beamsplitter  reflects  roughly  50%  (Fig.  1).  Slight  adjustments  of  a  camera's  location  can  be  used 
to  fine-tune  peak  intensity  siiKe  maximum  intensity  drops  as  the  laser  beam  diverges. 

The  integration/storage  capacity  of  CCD  arrays  combined  with  the  ability  to 
simultaneously  capture  a  video  frame  from  three  independent  cameras  allows  a  measurement  to 
be  made  from  a  single  pulse  of  laser  light,  once  the  waist  location  has  been  found.  The  Sony 
SSC-D5  cameras  are  locked  in  synch  and  run  continuously.  A  pulse  of  light  from  the  laser 
triggers  the  Data  Translation  CJuickCapture  video  frame  grabber  boards  to  capture  the  next 
frame  in  the  video  sequence  from  each  camera.  Those  frames,  which  each  of  the  three  boards 
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converts  into  an  image,  contain  the  transverse  laser  beam  profiles  at  the  cameras.  Each  image  is 
composed  of  640  x  480  pixels  using  256  grey  levels.  Laser  pulses  of  interest  will  in  general  be 
short  compared  to  the  single  frame  integration  time  of  the  CCD  cameras  (17  ms). 

The  effective  beam  widths  can  then  be  determined  in  software  in  any  of  several 
different  ways,  such  as  simulation  of  a  knife  edge  scan  (for  any  clip  level),  which  is  very  fast, 
or  direct  computation  of  variances,  which  is  more  theoretically  correct.  Fast  assembly  language 
routines  are  used  for  all  computations,  including  determination  of  principal  axes,  to  provide 
inunediate  feedback  on  beam  quality.  The  beam  quality  is  determined  from  the  product  of  the 
standard  deviation  of  the  transverse  beam  profile  at  the  beam  waist  gq,  and  the  standard 
deviation  Gz,  of  the  beam  profile  at  a  known  distance  away,  z,  in  the  form^ 

zX 

The  two  far-field  cameras  provide  two  measurements  of  the  value. 

It  would  also  be  possible  to  measure  beam  quality  by  finding  the  best-fit  parabola 
through  the  beam  widths  of  all  three  transverse  profiles.  However,  this  method  is  very 
sensitive  to  small  errors  in  the  size  of  the  beam  waist.  Thus,  it  is  much  better  to  find  the  beam 
waist  independently  and  calculate  the  best-fit  parabola  while  holding  the  beam  waist 
constant.  In  our  case,  this  would  amount  to  little  more  than  averaging  the  two  values  of  as 
obtained  by  the  two  cameras  as  outlined  above. 

We  have  tested  the  operation  of  this  system  using  single  pulses  from  a  HeNe  laser  with 
various  artificial  distortions  added,  and  are  now  beginning  a  series  of  measurements  on  other 
pulsed  laser  sources.  To  simulate  a  distorted  pulsed  laser  source,  a  CW  HeNe  laser  beam  was 
sent  through  an  acousto-optic  modulator  and  a  prism  pair.  The  AO  modulator  produces  62  ps 
square-wave  pulses,  and  the  prism  pair  elongates  and  slightly  distorts  the  beam  along  one  axis. 
Figure  2  shows  the  three  transverse  profiles  taken  during  a  single  pulse.  The  principal  axes  and 
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beam  widths  as  measured  by  the  computer  have  been  added  to  the  figure.  The  lengths  of  the 
ellipse  semiaxes  are  equal  to  the  laser  beam  image  variances.  The  principal  axes,  as  calculated 
by  the  computer,  differ  from  each  other  by  about  one  degree.  We  assume  that  this  is  due  to 
small  alignment  errors  and  not  to  the  laser  beam  being  nonorthogonal^.  The  fringes  seen  in  the 
transverse  profiles  do  not  greatly  effect  the  value  of  if  the  power  in  the  fringes  is  relatively 
small.  The  astigmatism  of  the  measured  laser  beam  was  negligable. 

The  distance  from  the  beam  waist  to  camera  B  in  this  case  was  zb  =  124  ±  1  mm  (Fig.  1). 
The  distance  from  the  beam  waist  to  camera  C  was  Zc  =  173  ±  1  mm.  A  single  pixel  measures  13.4 
^m  on  a  side.  At  the  beam  waist  the  pixels  of  camera  A  are  the  same,  but  a  40x  microscope 
objective  makes  the  effective  length  0.335  pm  on  a  side.  Repeated  measurements  on  different 
pulses  from  the  same  laser  give  variances  which  differ  by  roughly  one  pixel. 

The  measured  beam  widths  were: 

Camera  A  (waist):  OQx  =  75  pixels  (25.1  pm)  t^Oy  =  21  pixels  (7.04  pm) 

Camera  B:  Ocx  =  23  pixels  (308  pm)  Ocy  =  67  pixels  (898  pm) 

Camera  C:  Obx  =  32  pixels  (429  pm)  oby  =  97  pixels  (1300  pm) 

From  this  data,  using  the  previous  equation,  the  calculated  m2  values  are; 
can^a  A  (waist)  and  camera  B:  Mx^  =  1-24  My^  =  1.01 

camera  A  (waist)  and  camera  C:  Mx^  =  1.24  My^  =  1.05 

The  two  cameras  agree  to  within  about  5%  of  each  other. 

We  expect  to  extend  these  measurements  to  other  pulsed  laser  sources,  such  as  pulsed 
semiconductor  diode  lasers  and  CJ-switched  diode-pumped  lasers.  Given  the  interesting 
phenomena  already  found  when  examining  the  beam  quality  of  cw  laser  sources  for  different 
ranges  of  operation,3  we  expect  to  find  similarly  interesting  results  in  the  modal  and  beam- 
development  properties  of  pulsed  laser  sources. 
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Figure  Captions 

Fig.  1.  Laboratory  setup  for  measuring  Beam  Quality  of  Pulsed  lasers. 
Fig.  2.  Transverse  beam  profiles  of  a  single  laser  pulse. 
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ABSTRACT 

A  recent  reformulation  of  the  laser  field  equations,  using  the  true  modes  of  the  laser  cavity  and 
including  a  noise  polarization  term  to  account  for  spontaneous  emission,  shows  that  conventional 
understanding  of  amplifier  noise  figure,  oscillator  linewidth,  and  oscillator  build-up  time  needs  to  be 
modified  for  laser  amplifiers  or  oscillators  with  nonorthogonal  transverse  eigenmodes.  Gain-guided 
lasers  and  unstable-resonator  lasers  in  particular  have  nonorthogonal  transverse  modes,  and  hence 
exhibit  unconventional  noise  behavior.  One  important  consequence  of  this  theory  is  the  appearance 
of  an  excess  spontaneous  emission  factor  (the  “Petermann  factor”)  which  multiplies  the  well-known 
Schawlow-Townes  formtila  for  the  laser  linewidth.  This  factor  can  be  substantially  greater  than  imity 
for  unstable  resonators  with  large  magnification  or  FVesnel  number.  We  discuss  our  investigations  of 
the  consequences  of  the  nonorthogonality  of  the  transverse  eigenmodes,  especizilly  on  laser  linewidth 
in  tmstable-resonator  lasers.  We  believe  that  the  results  of  these  investigations  are  important 
where  the  temporal  coherence  properties  of  the  laser  are  of  interest,  namely  in  coherent  optical 
communications,  spectroscopy,  optical  wavelength  standzu’ds,  and  also  in  the  injection  seeding  of 
pulsed  high-power  unstable  resonator  lasers.  Extensive  calculations  of  the  excess  noise  factor  in  real 
hard-edged  unstable  resonators  have  been  carried  out  using  a  virtual-source  approach  to  calculate 
the  exact  unstable-resonator  eigenmodes.  An  experiment  is  also  being  carried  out  to  demonstrate 
the  linewidth  enhancement  in  a  small  diode-pumped  geometrically  unstable  laser  oscillator. 

1.  INTRODUCTION 


In  many  applications,  understeinding  laser  linewidth  is  an  import  Eint  consideration  for  example 
for  the  design  of  coherent  optical  communication  schemes,  for  spectroscopy  (when  investigating 
fundamental  properties  of  materials),  and  for  the  use  of  lasers  in  gravity  wave  detection  using 
interferometers.  The  fundamenteil  source  of  laser  linewidth  is  spontaneous  emission  noise  which  is 
the  unavoidable  process  accompanying  stimulated  emission.  It  is  this  contribution  to  laser  linewidth 
we  are  interested  in  as  opposed  to  external  noise  contributions  which  depend  on  the  environment 
in  which  the  ideal  laser  with  spontaneous  emission  noise  is  imbedded. 

Let  us  now  turn  our  attention  to  laser  resonators.  It  is  well  known  that  unstable  resonator 
lasers  have  several  advantages  over  stable  resonator  lasers,  in  paurticular  good  transverse  mode  dis¬ 
crimination  and  large  gain  volume.  Two  examples  of  practical  lasers  built  with  unstable  resonators 
are  CO2  lasers  and  some  Q-switched  Nd:YAG  lasers.  The  property  that  we  are  particulaxly  inter¬ 
ested  in  is  that  the  modes  of  unstable  resonators  are  nonorthogonal  eind,  as  we  shall  describe  in  this 
paper,  this  has  direct  bearing  on  the  fimdamental  linewidth  of  the  laser. 

Note  that  for  most  lasers,  the  environmental  or  external  noise  sources  typicailly  dominate  the 
contribution  due  to  spontaneous  emission.  An  important  exception  is  the  semiconductor  Iztser  whose 


1 


linewidth  is  dominated  by  spontaneous  emission  noise.  We  can  therefore  anticipate  that  by  building 
semiconductor  lasers  with  unstable  resonators  or,  if  more  generally  these  semiconductor  lasers  have 
nonorthogonal  eigenmodes,  they  will  exhibit  excess  spontaneous  emission  noise  properties  which 
bear  directly  on  their  fundamental  linewidth. 

Note  that  when  we  say  “quantum  linewidth”  or  “Schawlow- Townes  linewidth”  or  “fundamental 
linewidth”  of  a  laser  we  mean  the  spontaneous  emission  noise  linewidth. 

We  shall  first  review  the  orthogonality  properties  of  modes  in  lasers,  then  outline  the  theory  of 
excess  quantum  noise,  and  give  some  numerical  results  for  the  excess  quantmn  noise  for  two  types 
of  lasers,  zind  follow  by  a  discussion  of  experimental  issues  involved  in  trying  to  measure  the  excess 
quantum  noise. 


2.  MODES  IN  RESONATORS 

We  wish  to  give  some  details  concerning  resonator  modes  and  their  properties  in  order  to  avoid 
referring  the  reader  to  a  large  collection  of  papers  discussing  one  or  another  aspect  of  the  topic. 
(See  reference  1  for  background.) 

Consider  the  unstable  optical  resonator  in  Fig.  1.  We  define  a  reference  plane  in  the  resonator 
which,  for  convenience,  will  be  taken  as  the  plane  transverse  to  the  optic  axis  just  before  the  output 
coupling  mirror.  Note  that  in  the  paraxial  approximation,  the  output  mirror  is  considered  to  be  of 
negligible  width.  We  meein  by  “mode”  or  “eigenmode”  of  this  resonator  any  field  distribution  which 
will  reproduce  after  one  round  trip  through  the  resonator.  In  order  to  make  the  picture  clearer,  it  is 
useful  to  consider  the  following  lensguide  equivalent  where  the  resonator  is  unfolded  and  the  mirrors 
are  replaced  by  lenses.  Also,  coupling  to  the  outside  is  denoted  by  the  presence  of  an  absorbing  wall 
(representing  radiation  going  out  to  infinity)  as  shown  in  Fig.  2.  A  field  distribution  u  at  plane  Zi  v., 
an  eigeiunode  if  when  propagated  to  plane  Z2  (one  round-trip  through  the  resonator)  it  reproduces 
its  shape  both  in  amplitude  and  phase,  so  that 

vis,  Z2)  =  l  u{s,Zi) 

where  s  is  the  treinsverse  dimension.  Note  that  u  is  the  slowly  varying  part  of  a  field 

The  constant  7  is  the  eigenvalue  associated  with  the  eigenmode  u.  The  mode  u  obeys  the  paraxial 
wave  equation  with  appropriate  boundairy  conditions 

In  O’^der  to  understand  the  orthogonality  property  of  modes  it  is  more  instructive  to  consider 
the  integral  equation  obeyed  by  u.  This  is  nothing  else  but  the  Fresnel  diffraction  integral 

^2('S2)  =  /  ■A(S2,o>i)U](Si)dSi 

Ja 

where  A  is  the  aperture  area,  and  where  the  kernel  K  depends  on  the  optical  elements  between  the 
reference  planes  at  zj  ^lnd  22  and  the  integral  is  taken  over  tlie  aperture  in  the  equivalent  lensguide. 
We  write  Ui(si)  as  a  shorthand  for  ui(si,Zj),  and  likewise  for  V2-  For  simplicity  of  argument, 
we  will  consider  the  one- dimensional  case  where  s  =  x  and  S2  =  y.  This  choice  is  by  no  means 
restrictive,  but  it  makes  the  notation  easier  and  also,  for  the  two-dimensional  case,  it  suffices  to 


Fig.  2.  Lensguide  equivalent  of  the  resonator  shown  in  Fig.  1.  The  arrows  represent  the  to* ward 
propagating  eigenwave.  The  reference  plcines  1  and  2  are  located  just  before  the  absorbing  walls. 


view  X  as  the  vector  {x\^X2)  and  y  as  the  vector  (yi,y2)-  So  we  write  the  integral  equation  for  the 
mode  as 

v{y)  =  J  K{y,x)u(x)dx 

A  mode  of  the  resonator  is  the  self-reproducing  field  u  such  that 

v{y)  =  ju(y) 
or 

7u(y)  =  J  K{y,x)u{x)dx 

so  that  u  is  an  eigenmode  of  the  integral  operator  which  transforms  u(x)  into  J  K(y,x)u(x)dx 
with  its  corresponding  eigenvalue  7. 

Let  us  suppose  that  {u„},  {7n}  is  the  collection  of  eigenmodes  and  corresponding  eigenvalues. 
We  define  the  scalar  product  hy 

(Um,Un)  =  J  U^{x)Unix)dx 

where  the  integral  is  talcen  over  the  aperture.  In  many  cases  in  physics  the  eigenmodes  are  orthog¬ 
onal,  i.e. 

(I  n  =  m 

(Unn^n)  =  ^mn  =  1  ^  . 

10  n  ^  m 

but  in  the  case  of  unstable  resonators  this  is  not  the  case. 


{um,Un)  ^  Smn  for  unstable  resonators 

It  is  possible  however  to  define  a  set  of  modes  {um}  which  is  orthogonal  to  the  set  of  modes  {um} 

=  ^nm 


It  turns  out  that  these  modes  are  the  eigenmodes  of  the  integral  equation 


Im^m(y)  =  J  K*{X,  y)Um(l)  dx 


Note  that  the  arguments  of  K  are  switched  and  it  is  complex  conjugated.  The  field  Vm  is  called  an 
adjoint  mode  of  Um.  In  laser  resonators  another  definition  is  used:  we  define  the  mode  (j>m  as  the 
solution  to 

K(x,y)4)m{x)dx  (no  complex  conjugate) 

So  that  and  it  obeys  the  biorthogonality  property 


—  ^mn 


The  reason  for  using  this  definition  is  that  <l>m  can  be  interpreted  as  the  eigenmode  of  the  resonator 
propagating  in  the  opposite  direction  (with  the  optical  axis  reversed). 

It  turns  out  that  for  optical  resonators  in  the  peiraxial  approximation  we  have 

<f>mix)  =  /9m  Um(a:)e“’“** 

where  the  constant  a  is  given  for  a  given  resonator.  If  the  reference  plane  is  a  symmetry  plane  of 
the  resonator  we  have  a  =  0,  and  then  the  biorthogonality  relation  becomes 

J  Um(x)un{x)  dx  =  0  m^n  (no  complex  conjugation) 

It  is  necessary  also  to  define  a  normalization  of  the  modes.  We  will  choose 

/  ltim(a:)!^  di  =  1 
J  — OO 

and 

/  <f)m{x)u^{x)dx  =  1 

Ja 

But  from  the  previous  equation  we  have  that 

/  /9mUm(a;)e“*“**Um(x)dx  =  1 

Ja 


so 


/9m  = 


Ja  ult(x)e-‘°^"  dx 


It  is  easy  to  show  that 


/  lum(x)|^di  =  l7ml  ^  f  lum(a:)|^dx 

9-00  Ja 

so  we  have  fin2illy  a  set  of  eigenmodes  {um}  with  corresponding  eigenvedues  {7m}  such  that 

7mWm(y)  =  J  K{y,x)Um{x)dx 

with  the  normalization 

/  |«m(a:)pdx  =  |7m|^ 

Ja 


and  (^m  such  that 


<l>m{x)  = 


Um{x)c 
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The  quantity  we  axe  ultimately  interested  in  is  the  excess  noise  factor  Kjnm  defined  by 


Ja 


(x)<f>mix)dx 


And  from  the  above  relations  we  have 


r-  Ia  u^(x)um(x)dx 

•^mm  —  .  ^  ^ 

IL  dx\ 


which  at  a  symmetry  plane  simplifies  to 


A'„„  =  L  l"m(x)pA 


3.  EXCESS  NOISE  FACTOR  THEORY 

We  shall  describe  in  words  the  derivation  leading  to  the  excess  noise  factor.  For  the  full  theory 
we  refer  the  reader  to  reference  2. 

The  gist  of  the  theory  is  to  go  from  Maxwell’s  equations  for  the  electric  field  to  the  equation 
for  the  photon  number  in  a  single  lasing  mode  by  an  expansion  of  the  field  in  the  true  resonator 
eigenmodes  (longitudinal  and  transverse).  Conventional  theories  assume  am  orthogonal  mode  de¬ 
composition  using  the  modes  of  a  closed  cavity  with  no  losses.  The  laser  resonator  is  inherently 
an  open  resonator  and  cannot  be  correctly  approximated  by  a  closed  resonator  especially  for  lasers 
with  diffractive  output  coupling  such  as  provided  by  imstable  resonators. 

The  central  result  we  are  interested  in  is  the  equation  of  motion  for  the  photon  number  in  a 
single  Izising  mode 

dn 

—  =  KN2{n  -h  K)  —  kNiU  —  7cn 
dt 

We  recognise  the  familiar  stimulated  emission  term  KN2'n  (k  is  a  coupling  constant),  the  stimulated 
absorption  —kNiti,  emd  the  photon  loss  7cn.  The  spontaneous  emission  is  given  by  KN2K  which 
is  different  from  the  conventional  result  KN2  {K  replaced  by  1).  The  factor  K  depends  directly  on 
the  orthogonality  property  of  the  modes.  It  is  composed  of  a  “longitudinal”  part 


and  a  “transverse”  pcirt 


'  -  W  1  -  i7or  \ 
M  Vln(l/|7oP)/ 

AV  =  J  (f>l{x)<i)„{x)  dx 


where  <f>o  is  the  backward-propagating  mode  associated  with  the  lasing  mode  Uo- 

The  factor  A'/,  is  close  to  unity  except  for  resonators  with  very  large  output  coupling.  It  is 
based  on  the  nonorthogonadity  of  the  longitudinal  modes  of  the  laser.  The  factor  A'r  is  close  to  unity 


for  stable  resonator  lasers  and  can  be  much  larger  than  unity  (as  we  shall  see  later)  for  unstable 
resonator  lasers  which  have  strongly  nonorthogonal  eigenmodes. 


1.  QUANTUM  LINEWIDTH 


The  Schawlow- Townes  linewidth  or  quantum  linewidth  of  the  laser  is  determined  by  the  amoimt 
of  spontaneous  emission  into  the  lasing  mode.  If  the  laser  has  nonorthogonal  eigenmodes  we  conclude 
that  the  amount  of  spontaneous  emission  in  the  lasing  mode  is  larger  by  a  factor  K  than  predicted 
by  a  conventional  laser  theory  based  on  an  orthogonal  mode  expansion.  We  therefore  expect  that 
the  qu2mtum  linewidth  of  the  laser  will  be  increased  by  the  same  factor.  The  expression  for  the 
laser  linewidth  is  then 


A  I//, 


2-KhiiJ 

P 


■K 


The  first  two  terms  represent  the  conventional  laser  linewidth  depending  on  the  cold  cavity  linewidth 
Ai/c  and  the  output  power  P.  We  should  note  that  in  general,  two  more  factors  may  appear  in 
the  above  expression,  namely  the  incomplete  inversion  factor  N2/{N2  —  Ni)  (important  for  3-level 
lasers)  euid  the  linewidth  enhancement  factor  (1  -f-a^)  for  lasers  exhibiting  AM  to  FM  noise  coupling 
such  as  semiconductor  lasers.^  The  excess  spontaneous  emission  factor  K  is  then  the  product  of  the 
longitudinal  factor  and  the  transverse  factor  Kt- 

The  longitudinal  {Ki)  factor  has  been  investigated  experimentally  by  Hamel  and  Woerdman,^ 
and  is  important  for  lasers  with  very  large  output  couplings.  We  are  particularly  interested  in  the 
transverse  (Kt)  factor  because  it  is  greater  than  unity  for  imstable  resonator  Izisers. 


5.  EXCESS  NOISE  FACTOR:  NUMERICAL  PREDICTIONS 

5.1.  Variable  reflectivity  mirror  laser 

We  have  investigated  the  case  of  variable  reflectivity  mirror  (VRM)  lasers  partly  because  such 
lasers  are  increasingly  popular.  Also,  in  the  ceise  of  a  gaussian  reflectivity  profile,  the  eigenmodes 
and  eigenvalues  can  be  computed  analytically.  Detailed  computations  can  be  found  in  reference 
5,  but  we  sh^lll  summarize  the  essential  features  here.  Figure  3  shows  a  plot  of  the  excess  noise 
factor  of  the  lowest  loss  mode  in  a  VRM  resonator  laser  as  a  function  of  the  geometrical  stability 
parameter  m  (the  half-trace  of  the  round-trip  ABCD  matrix  representing  the  resonator).  In  the 
stable  resonator  region  (|m|  <  1)  the  excess  noise  factor  is  essentially  unity  so  that  there  is  no 
excess  noise.  For  a  flat-flat  resonator  (m  =  1),  the  excess  noise  factor  is  equal  to  2,  and  in  the 
unstable  resonator  region  the  excess  noise  factor  increeises  very  rapidly  and  reaches  a  value  of  100 
at  a  magnification  M  (M  =  m  -f  —  1)  of  1.6  for  this  particular  case.  We  conclude  that  VRM 
unstable  resonator  lasers  can  exhibit  very  large  excess  noise  factors. 

5.2.  Hard-edged  unstable  resonator  laser 

In  this  case,  the  eigenmode  problem  cannot  be  solved  analytically,  but  we  have  computed 
the  excess  noise  factor  using  an  efficient  algorithm  called  the  virtual  source  method  which  allows 
generation  of  plots  of  the  excess  noise  factor  without  having  to  rely  on  standard  algorithms  (the 
Fox- Li  method  or  Prony  method)  which  would  be  much  more  time-consuming. 
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Fig.  3.  Power  loss  1  —  I7P  in  percent  per  round  trip  and  excess  noise  factor  K  plotted  versus 
geometrical  magnification  for  the  lowest-order  Hermite-gaussian  mode  in  a  VRM  cavity  with  a  ratio 
of  gaussian  aperture  to  confocal  spot  size  of  Wgalwcf  =  5  or  gaussian  Fresnel  number  Nga  «  8. 
The  horizontal  scale  variable  in  the  half-trace  parameter  m  for  0  <  m  <  1,  and  the  geometrical 
magnification  Af  =  m  -f  x/m^  —  1  for  m  >  1. 


EFFECTIVE  FRESNEL  NUMBER 


Fig.  4.  Power  loss  and  excess  noise  factor  versus  effective  Fresnel  number  ter  a  hard-edged  strip 
unstable  resonator  with  a  magnification  of  4. 


We  plot  here  power  losses  and  the  excess  noise  factor  as  a  function  of  the  effective  Fresnel 
number  in  the  resonator  (i.e.,  as  a  function  of  output  coupling  mirror  diameter,  all  other  parameters 
being  kept  constant).  Figure  4  shows  how  the  excess  noise  factor  is  “resonant”  with  the  Fresnel 
number  in  the  same  way  as  are  the  losses.  This  fact  is  typical  of  hard-edged  mirror  resonators. 
Again,  we  see  that  the  excess  noise  factor  can  reach  very  large  values. 

6.  EXPERIMENTAL  CONSIDERATIONS 

We  are  currently  setting  up  an  experiment  to  measure  the  excess  noise  factor  using  a  diode- 
pumped  monolithic  Nd;YAG  laser  system  (Fig.  5).  The  rod  has  a  flat  mirror  on  one  end  and  a 
negative  curvatiire  nurror  on  the  other.  It  is  thus  inherently  geometrically  unstable.  The  rod  is 
end-pumped  by  a  single  mode  high-power  diode  laser  which  induces  a  lasing  mode  in  the  rod  by  the 
mechanism  of  gain-guiding.  Although  the  theory  presented  above  makes  the  assumption  that  the 
gain  is  uniform  throughout  the  laser,  we  anticipate  an  effect  on  the  excess  of  noise  from  the  mirror 
curvature. 

For  the  noise  measurement  we  are  setting  up  a  Pound-Drever-type  noise  measurement  scheme® 
where  a  Fabry-Perot  cavity  is  locked  to  the  laser  and  frequency  noise  is  observed  at  frequencies 
beyond  the  locking  bandwidth,  i.e.,  where  the  laser  is  free-running  and  where  external  noise  sources 
have  become  negligible,  leaving  the  spontaneous  emission  noise. 

We  have  observed  lasing  action  in  our  diode-pumped  rod  system  at  a  rod  temperature  of  80K 
when  the  gain  is  high  enough  to  go  above  threshold  and  to  avoid  multi-longitudinal  mode  oscillation 
by  narrowing  of  the  gain  curve  (hence  avoiding  spatial-hole  burning). 

We  wish  to  stress  that  this  laser  scheme  is  interesting  in  itself  because  it  has  good  transverse 
mode  discrimination  in  a  large  gain  volume,  so  that  it  is  not  limited  by  the  problem  of  stable  rod 
designs  where  aJl  the  pump  power  has  to  be  mode  matched  into  the  fundamentEil  mode  of  the  rod 
Iciser.  Figure  6  shows  the  near-field  beam  emanating  from  the  rod. 

7.  CONCLUSIONS 

Unstable  resonator  Ijisers  can  exhibit  large  deviations  from  standard  predictions  for  the  fun¬ 
damental  leiser  linewidth  due  to  eigenmode  nonorthogonality.  Indeed,  the  standard  expression  for 
the  laser  linewidth  due  to  spontaneous  emission  noise  needs  to  be  multiplied  by  a  so-called  “excess 
spontaneous  factor.”  We  have  performed  analytic  and  numericaJ  computations  for  the  excess  noise 
factor,  and  conclude  that  it  can  be  very  much  larger  than  unity  for  unstable  resonators.  Note  that 
the  excess  noise  factor  theory  was  first  studied  in  the  context  of  gain-guided  semiconductor  lasers 
by  Petermann.^  Therefore,  it  is  also  called  the  Petermann  factor. 

We  are  currently  investigating  the  excess  noise  factor  experimentally,  and  have  demonstrated 
what  is  to  our  knowledge  the  first  diode-pumped  monolithic  unstable  resonator  Nd:YAG  laser 
system. 


Pump  beam 


NdiYAG  laser  rod  Output  laser  beam 


Fig.  5.  Optical  set-up  of  the  diode-pumped  monolithic  unstable  resonator  Nd:YAG  rod  laser. 
The  rod  is  3  mm  long  with  a  flat  mirror  at  one  end  (HR  at  1064  nm  and  HT  at  809  nm)  and  a  —10 
cm  radius  mirror  at  the  other  end  (95%  reflector  at  1064  nm).  The  material  is  1%  doping  Nd:YAG. 


Fig.  6.  Near  field  of  the  bearr  emanating  from  the  rod.  On  the  left  the  laser  is  below  threshold 
and  on  the  right  the  laser  is  above  threshold.  Notice  a  change  in  vertical  scale  by  a  factor  of  1000 
and  a  change  in  horizontal  scale  by  a  factor  of  3/4.  The  ring  of  light  in  the  left  pattern  is  the  glow 
from  the  edge  of  the  laser  rod. 
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